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1  Executive  Summary 

This  report  presents  the  results  of  the  work  performed  with  unflawed  notched  NACA  specimens 
used  in  connection  with  Contract  No.  N68335-10-C-0428  -  Option  1,  Task  1.3. 

The  integrated  software  system  and  computational  framework  developed  by  ESRD  under  the 
Navy  SBIR  2008.2  Phase  II  project  -  Topic  N08-131:  Innovative  Approaches  to  Flaw-Tolerant 
Design  and  Certification  of  Airframe  Components  is  a  new  technology  ideally  suited  to  model 
and  analyze  structural  damage  that  exceeds  existing  repair  limits;  to  rank  in-service  damage 
based  on  fatigue  criticality,  not  just  damage  size;  to  enable  the  quantitative  assessment  of  risks  to 
flight  safety  from  using  condition-based  maintenance  (CBM)  to  delay  or  eliminate  fleet 
inspection  requirements;  to  meet  CBM  need  to  move  flaw-damaged  component  removals  from 
unscheduled  to  scheduled  maintenance  action  by  enabling  prediction  of  remaining  life. 

One  of  the  unique  aspects  of  the  approach  is  the  development  of  a  family  of  damage 
accumulation  models  for  ranking  in-service  damage  based  on  fatigue  criticality.  The  proposed 
stress-based  damage  accumulation  models,  using  the  integral  average  of  a  function  of  stress  over 
a  solution-dependent  volume,  were  tested  against  available  experimental  data  available  in  the 
literature  for  notched  fatigue  specimens  and  with  flaw-seeded  fatigue  data  provided  by  the  US 
Navy.  The  characterizing  parameter  of  the  models  of  damage  accumulation  if)  was  correlated 
with  the  highly  stressed  volume  ( HSV)  of  material  around  the  flaw  to  allow  the  generalization  of 
the  calibration  work  performed  for  flaw  seeded  specimens  to  components  with  flaws  of  various 
sizes  and  shapes. 

This  report  presents  the  work  performed  with  unflawed  notched  NACA  specimen  and  provides 
step  by  step  procedures  to  obtain  the  reference  S-N  curve  from  unnotched  specimens  of  24S-T3, 
75S-T6  and  SAE  4130,  the  calibration  procedure  to  obtain  the  J3-HSV  curves  for  each  material, 
and  the  prediction  of  fatigue  life  based  on  the  calibration  curve,  including  the  measure  of  the 
difference  between  model  prediction  and  the  experimental  results.  The  report  includes  all  data 
extracted  from  the  relevant  NACA  reports,  which  is  available  in  the  SimGov  Data  Management 
Portal  https://ftd.esrd.com. 
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2  Task  6  -  Report  on  NACA  Data 

2.1  Task  description 

The  following  description  corresponds  to  Task  6  of  the  Option  3  Period  of  the  contract: 

•  Prepare  a  detailed  report  of  the  work  performed  with  unflawed  notched  NACA  specimens 
used  in  connection  with  Contract  No.  N68335-10-C-0428  -  Option  1,  Task  1.3.  Task  will 
document  the  step  by  step  procedure  to  obtain  the  reference  S-N  curve  from  unnotched 
specimens  of  24S-T3,  75S-T6  and  SAE  4130,  the  calibration  procedure  to  obtain  the  [3- 
HSV  curves  for  each  material,  and  the  prediction  of  fatigue  life  based  on  the  calibration 
curve,  including  the  measure  of  the  difference  between  model  prediction  and  the 
experimental  results.  The  report  shall  include  all  data  extracted  from  the  relevant  NACA 
reports  available  in  the  SimGov  Data  Management  Portal  https ://ftd.  esrd.  corn 

2.2  Summary  of  NACA  data 

The  NACA  data  consists  of  a  collection  of  technical  notes  containing  axial-load  fatigue  test 
results  for  unnotched  and  notched  sheet  specimens  of  24S-T3,  75S-T6  aluminum  alloys  and  SAE 
4130  steel,  as  well  as  all  the  tabulated  data  in  spreadsheets  deployed  in  the  SimGov  Data 
Management  portal  (see  Appendix  4  for  details).  The  references  for  the  technical  notes  are  as 
follows: 

•  Grover  HJ,  Bishop  SM  and  Jackson  LR,  Fatigue  Strength  of  Aircraft  Materials.  Axial-Load 
Fatigue  Tests  on  Unnotched  Sheet  Specimens  of24S-T3  and  7 5 S-T6  Aluminum  Alloys  and  of 
SAE  4230  Steel.  NACA  Technical  Note  2324,  March  1951. 

•  Grover  HJ,  Bishop  SM  and  Jackson  LR,  Fatigue  Strength  of  Aircraft  Materials.  Axial-Load 
Fatigue  Tests  on  Notched  Sheet  Specimens  of  24S-T3  and  75S-T6  Aluminum  Alloys  and  of 
SAE  4230  Stainless  Steel  with  Stress  Concentrations  Factors  of  2.0  and  4.0.  NACA 
Technical  Note  2389,  June  1951. 

•  Grover  HJ,  Bishop  SM  and  Jackson  LR,  Fatigue  Strength  of  Aircraft  Materials.  Axial-Load 
Fatigue  Tests  on  Notched  Sheet  Specimens  of  24S-T3  and  75S-T6  Aluminum  Alloys  and  of 
SAE  4230  Stainless  Steel  with  Stress  Concentrations  Factors  of  5.0.  NACA  Technical  Note 
2390,  June  1951. 

•  Grover  HJ,  Bishop  SM  and  Jackson  LR,  Fatigue  Strength  of  Aircraft  Materials.  Axial-Load 
Fatigue  Tests  on  Notched  Sheet  Specimens  of  24S-T3  and  75S-T6  Aluminum  Alloys  and  of 
SAE  4230  Stainless  Steel  with  Stress  Concentrations  Factors  of  1.5.  NACA  Technical  Note 
2639,  February  1952. 
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•  Grover  HJ,  Bishop  SM  and  Jackson  LR,  Fatigue  Strength  of  Aircraft  Materials.  Axial-Load 
Fatigue  Tests  on  Notched  Sheet  Specimens  of  2024-T3  and  7075-T6  Aluminum  Alloys  and  of 
SAE  4230  Stainless  Steel  with  notched  radii  of  0.004  and  0.070  inch.  NASA  Technical  Note 
D-lll,  September  1959. 

The  material  properties  and  thicknesses  of  the  NACA  specimens  are  shown  in  Table  1,  while  the 
details  of  the  test  records  are  summarized  in  Table  2  (values  in  brackets  in  the  last  3  columns  are 
the  number  of  specimens  that  failed  during  the  fatigue  test). 

Table  1:  Material  Properties  of  NACA  specimens 


Material 

Thickness 

(in) 

Grain 

Direction 

Average  Tensile  Propoerties 

Average  Compressive  Properties 

Elongation  in 
2-in  (%) 

Yield  Strength 
0.2%  Offset 

(psi) 

Ultimate 
Strength  (psi) 

Yield  Strength 

(psi) 

Modulus  of 
Elasticity  (psi) 

24S-T3 

0.090 

with 

18.20 

54,000 

73,000 

44,500 

10,650,000 

cross 

18.30 

50,000 

71,000 

50,500 

10,450,000 

75S-T6 

0.090 

with 

11.40 

76,000 

82,500 

74,000 

10,450,000 

cross 

11.00 

75,000 

82,500 

78,500 

10,550,000 

SAE  4130 

0.075 

with 

14.25 

98,500 

117,000 

86,000 

30,400,000 

cross 

12.50 

101,000 

120,000 

97,000 

31,300,000 

Table  2:  Summary  of  the  928  available  fatigue  test  records 
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2.3  Reference  S-N  Curve:  Implementation  of  Random  Fatigue  Models 

Mathematical  models  devised  for  the  prediction  of  damage  accumulation  in  high  cycle  fatigue 
consist  of  a  problem  of  linear  elasticity  plus  two  sub-models:  a  predictor  V  and  a  statistical 
model  S.  A  predictor  is  a  functional  of  the  solution  of  the  problem  of  linear  elasticity  (or,  more 
generally,  small-strain  solid  mechanics)  intended  for  the  generalization  of  the  results  of  uniaxial 
fatigue  data,  also  known  as  S-N  data,  to  triaxial  stress  conditions.  Details  are  available  in  [1], 

Statistical  models  are  formulated  for  the  quantification  and  generalization  of  the  statistical 
dispersion  of  S-N  data.  The  information  collected  from  constant  cycle  fatigue  tests  consists  of  (a) 
the  maximum  test  stress  Smax,  (b)  the  cycle  ratio  R,  which  is  defined  as  the  ratio  of  the  minimum 
test  stress  Smm  to  the  maximum  test  stress  Smax,  (c)  the  number  of  cycles  at  which  failure 
occurred  or  the  test  was  stopped  Nf  ,  (d)  an  indication  whether  the  test  was  stopped  prior  to 
failure  (runout)  and  (e)  notes  indicating  unusual  events  such  as  buckling,  failure  at  the  grips  or 
failure  outside  of  the  test  section  that  may  have  occurred.  Typical  statistical  models  used  in  the 
analysis  of  fatigue  data  and  their  properties  are  described  in  [2],  The  implementation  of  two 
random  fatigue  limit  models  is  described  in  this  document.  These  models  are  based  on  the  papers 
of  Pascual  and  Meeker  [3],  [4],  Additional  discussion  and  analysis  are  available  in  [2], 

2.3.1  Formulation  of  statistical  model  S-\ 


We  define  the  equivalent  stress  aeq  as  [5]: 


_  rt 

^  eq  ^  max 


1  -R 

V~ 2~  j 


It  is  assumed  that  logio  Nf  is  normally  distributed: 

log l0Nf  ~  N (ju  ,s  ) 


Eq.  1 


Eq.  2 


The  mean  ju  is  defined  by 

//  =  A  -  A  log io(c%  -  As),  0  <  A3  <  creq  Eq.  3 

where  A3  is  the  fatigue  limit.  In  random  fatigue  limit  models  A3  is  a  random  variable.  We  assume 
that  logio  A3  is  normally  distributed: 

log10  A3  Eq.  4 


We  assume  that  the  standard  deviation  s  is  constant.  Therefore  the  statistical  model  is 
characterized  by  five  parameters  9=  {A\  A2  s  ju/Sf}.  Another  statistical  model  labeled  S2,  which 
has  six  parameters,  is  described  in  Appendix  1. 
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We  introduce  the  notation  w  =  logio  Nf,  v  =  logio  A3  and  invoke  two  fundamental  theorems  in 
statistics;  the  product  rule  and  marginalization.  The  product  rule  applied  to  the  random  fatigue 
limit  model  is: 

Pr  (vv,  v  I  aeq  )  =  Pr  (vv  I  <?eq ,  v)  Pr  (v  I  aeq  )  Eq.  5 

where  Pr  stands  for  the  probability  density  function.  The  equation  of  marginalization  is 

Pr  ( vv  I  aeq  )  =  |_J'°  Pr  (w,  v  I  <Jeq )  dv.  Eq.  6 

Therefore  the  probability  of  vv  given  aeq  is: 

Pr  (w  I  crcY/ )  =  J  J10  Pr(vv  I  <Jeq, v) Pr  (v  I  creq ) dv.  Eq.  7 

According  to  Eq.  (2)  w  has  normal  distribution  with  mean  ju  and  standard  deviation  5.  Therefore 
the  first  term  of  the  integrand  in  Eq.  (7)  is  the  probability  density  function: 


F(w,(jeq,v)  =  Pr(w  I  <Jeq,v)  = 


exp(-(  w  -  /Mcr^  ,10v  ))2  /(2^2 )) 

S'Jlx 


Eq.  8 


The  second  term  of  the  integrand  in  Eq.  (7)  follows  from  Eq.  (4).  By  hypothesis  the  probability 
density  of  v  is: 


/ (v)  =  Pr(v  I  ueq) 


exp(-(v  -  /Jf)2  /  2s) ) 


s  jr  •%/  ~2n 


Eq.  9 


Therefore  the  marginal  probability  density  function  of  vv,  given  aeq,  is: 

Fm  (w^eq)  =  Pr(w  I  <Jeq)  =  |  F  ( w,10'  )  f(v)dv. 
The  marginal  cumulative  distribution  function  (CDF)  of  w,  given  aeq,  is: 


Eq.  10 


pw  \  pit 

®M(W’°eq)  =  \_r  FM  (t,<Teq)dt  =  -J 


^OglO  °eq 

f  f 

1  +  erf 

J— oo 

V  v 

W-MiVeg’  1QV) 


J) 


f{v)dv.  Eq.  11 


Given  a  statistical  model  characterized  by  a  set  of  parameters  0,  and  a  set  of  independent 
observations  (vv„  cr"1),  (i  =  1,  2, .  .  . ,  m)  the  likelihood  function  is 


Eq.  12 


i= 1 


where 
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[  0  if  the  test  resulted  in  failure 

[l  if  the  test  was  stopped  before  failure 


Eq.  13 


The  maximum  likelihood  estimate  of  0,  denoted  by  6 ,  maximizes  L(0),  or  equivalently,  the  log 
likelihood  function  LL(0)  defined  by 


LL(0)  =  |](l-^)ln(FM  (»„<'))+  f>,  ln(l-4>„  (»„<>)). 

i= 1  i= 1 


Eq.  14 


2.3.2  Theoretical  design  curves 

The  theoretical  design  curve  <4(aeq)  is  the  set  of  points  (login  N,  ocq)  that  satisfy  the  condition 

ProbCA^  <  N  I  <yeq,6k)  =  p  Eq.  15 

where  k  identifies  the  statistical  model  and  p  represents  the  acceptable  probability  of  failure 
which  is  typically  a  small  number,  such  as  p  =  10  5.  Theoretical  design  curves  are  obtained 
through  application  of  statistical  models  well  beyond  their  range  of  calibration.  Furthermore, 
theoretical  design  curves  strongly  depend  on  the  choice  of  the  statistical  model.  Therefore 
considerable  aleatory  and  epistemic  uncertainties  exist.  Because  events  of  such  low  probability 
occur  very  rarely,  validation  of  statistical  models  with  respect  to  design  curves  is  not  feasible. 

The  setting  of  factors  of  safety,  which  are  to  account  for  aleatory  and  epistemic  uncertainties,  is 
burdened  by  the  same  uncertainties.  Therefore  subjective  judgment  in  setting  factors  of  safety 
cannot  be  avoided.  For  further  discussion  we  refer  to  [6]. 

2.3.3  Implementation  of  the  statistical  model 

The  implementation  is  based  on  Matlab  R2016a  (9.0.0.341360)  64  bit  (win64).  The  main  routine 
for  Statistical  Model  Si  is  a  Matlab  function  called  MMPDSmodellRFL(mat).m.  The  argument 
refers  to  the  material  type;  mat  =  1,  2,  3  corresponds  to  24S-T3,  75S-T6  aluminum  alloys  and 
SAE  4130  steel  respectively.  This  function  performs  the  following  operations: 

1.  Reads  the  material  data  identified  by  the  argument  mat.  The  data  are  in  the  file  NACA 
Unnotched  Specimen  4-13-16.xlsx. 

2.  Calls  the  Matlab  function  mle  (maximum  likelihood  estimate)  which  returns  parameters 
for  user-specified  statistical  models.  It  accounts  for  censored  data  (runouts).  The  PDF  and 
the  CFD  functions  are  called  custpdfRFLl.m  and  custcdfRFLl.m  respectively. 
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The  integrals  in  equations  (10)  and  (11)  have  to  be  evaluated  numerically.  In  the  test 
implementation  the  Matlab  function  integral  is  used  for  this  purpose  with  tolerance 
settings  'AbsTol '  =  10“12,  'RelTol '  =  10“12. 

The  maximum  likelihood  estimate  involves  the  solution  of  a  highly  nonlinear  problem  for 
the  parameters  of  the  statistical  model.  Typically  there  are  multiple  local  maxima.  The 
maximum  is  sought  in  the  neighborhood  of  user-provided  seed  values  of  the  parameters. 

A 

One  can  never  be  certain  that  the  global  maximum  of  LL(0),  denoted  by  LL(  6 ),  was 
found.  It  is  advisable  therefore  to  try  several  seed  values  and  select  the  set  of  parameters 
that  corresponds  to  the  maximum  value  of  the  log  likelihood  function  given  by  Eq.  (14). 
Estimation  of  seed  values  is  discussed  separately. 

3.  Computes  the  log  likelihood  function  LL{  6 )  and  the  Akaike  information  criterion  (AIC). 

4.  Plots  the  S-N  data  and  the  function  p  defined  by  Eq.  (3).  It  also  plots  the  theoretical 
design  curve  at  probability  level  p  =  pr.  For  the  definition  of  p  see  Eq.  (15).  The  design 
curve  parameters  A\,  A2,  A3  are  returned  in  vector  b.  The  pertinent  commands  are  as 
follows: 

%  Design  curve  at  probability  pr 
pr  =  0.001; 

[NAx,SA,G]  =  DesignCurvelRFL (mat ) ; 

Fnorm2  =  contour (NAx, SA, G,  [pr  pr]  ,  '  r' )  ; 
b  =  DesignCurveFit 1 (Fnorm2 )  ; 

An  example  is  shown  in  Figure  1 . 


103  104  105  106  107  108 


N 

Figure  1:  S-N  data  for  24S-T3  aluminum  alloy.  The  mean  corresponding  to  Model  Si  is  labeled  //i(aeq),  the 
theoretical  design  curve  corresponding  top  =  0.001  is  labeled  ^i(aeq).  Runouts  are  represented  by  the  open 

triangular  markers. 


Page  12 


ESRD,  Inc. 


SBIR  II  Option  3  -  Task  6 


N68335-1 6-C-0088 


Theoretical  allowable  value 


Given  a  theoretical  design  curve  and  the  number  of  design  cycles  N,  it  is  possible  to  compute  the 
theoretical  allowable  value  from 


k,L-rfr‘(«)  Ei“ 

For  example,  referring  to  Figure  1,  for  N  =  106  we  have  (aeq)aii  =  19.4  ksi. 

Tolerance  settings 

The  computations  were  performed  with  the  following  tolerance  settings: 

Display :  'off' 

MaxFunEvals:  1000000 
Maxlter:  1000000 
TolBnd:  1 . 0000e-06 
TolFun :  1 . 0000e-06 
TolX :  1 . 0000e-06 
GradOb j :  'on' 

DerivStep:  1 
FunValCheck:  'on' 

Seed  parameters 

The  maximum  likelihood  estimate  requires  seed  parameters  which  are  used  by  the  Matlab 
function  mle  as  the  starting  parameter  values  in  the  iterative  process  for  estimating  6. 

The  user  provides  an  upper  and  lower  bound  for  A3,  denoted  by  U  and  L  respectively.  These  are 
the  bounds  of  the  S-N  data  in  the  high  cycle  range.  The  interval  U  -  L  is  divided  into  m  sub¬ 
intervals  (for  example,  m  =  100). 

Letting 


A30  =  L  +  (i -  Y)(U  - L) / m,  i  =  1, 2,...,  m  +  l  Eq.  17 

we  find  A,0  and  Af  by  linear  regression  to  fit  the  data  to 

log10  Nj  =  A)0  -  Af  log10  (<reg  -  Af  )  Eq.  18 

where  j  identifies  all  data  points  (a^1 ,  Nj )  that  satisfy  the  criterion 

cA;1  -  A)'1  >  0.  Eq.  19 


Runouts  are  treated  in  the  same  way  as  other  points.  We  denote  the  number  of  data  points  by  n„ 
compute  the  variance  of  each  i. 


Eq.  20 
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and  identify  the  index  i  that  corresponds  to  the  smallest  variance.  We  denote  that  value  by  I.  To 
estimate  the  seed  value  corresponding  to  the  standard  deviation  for  model  Si,  we  compute  the 
sample  standard  deviation: 


S  = 


v «/ 


Li>g„  N,  -  (4"  - -  A!'1  logI0(<>  -A^f 


\l/2 


Eq.  21 


The  computations  are  performed  using  the  Matlab  function  SeedValuesLSl.m.  The  seed  values 
for  Model  Si  are: 

=  {A">  A"  s  log,„(A«>)  0.05}  Eq.  22 

where  the  number  0.05,  the  seed  value  for  Sf ,  was  determined  by  trial  and  error. 


Note:  All  the  matlab  functions  and  input  data  for  both  statistical  models  are  contained  in  the 
files:  Model  SI  Archived.zip  and  Model  S2  Archived.zip  (see  Appendix  4  for  details). 


2.3.4  Computed  model  parameters 


The  computed  statistical  parameters,  the  parameters  that  define  the  theoretical  design  curve,  the 
input  parameters  L,  U,  Sf  and  the  corresponding  output  parameters  consisting  of  the  log 

likelihood  estimate  LL(  0 )  and  the  Akaike  information  criterion  (AIC)  are  displayed  in  Figure  2. 


NACA  materials:  Statitical  Model  SI  | 

24S-T3 

75S-T6 

SAE  4130 

Parameters  for  the  design  curves  P=l/1000 

A1 

7.193 

6.336 

7.228 

24S-T3 

75S-T6 

SAE  4130 

A2 

1.992 

1.531 

1.764 

A1 

7.271 

6.411 

7.581 

A3 

22.06 

26.79 

54.09 

A2 

2.210 

1.737 

2.318 

s 

0.125 

0.115 

0.249 

A3 

15.615  1 

19.694 

40.602 

s_f 

0.0488 

0.0432 

0.0389 

Seed  and  output  parameters  | 

24S-T3 

75S-T6 

SAE  4130 

U 

35 

35 

75 

L 

15 

15 

45 

s_f 

0.05 

0.05 

0.05 

LL 

-576.73 

-881.37 

-477.71 

AIC 

1163.5 

1772.7 

965.4 

Parameters  for  the  design  curves  P=l/100000 1 

24S-T3 

75S-T6 

SAE  4130 

A1 

7.252 

6.393 

7.621 

A2 

2.271 

1.796 

2.480 

A3 

13.767  1 

17.629 

36.520 

Figure  2:  Unnotched  NACA  data.  Computed  and  input  parameters  for  statistical  model  1. 
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The  parameters  of  the  theoretical  design  curve  are  computed  by  the  function 
DesignCurveMainRFL(mat,p )  where  mat  identifies  the  material  and  p  is  the  probability  in  Eq. 
(15).  The  AIC  is  a  measure  of  the  goodness  of  fit  between  the  statistical  model  and  the  data.  The 
data  corresponding  to  AIC  was  computed  using  the  method  used  in  [2]  which  differs  by  scaling 
from  the  method  used  in  [4], 

The  goal  is  to  minimize  this  value,  given  the  data.  Since  the  statistical  model  depends  on  the 
values  of  the  seed  parameters,  the  seed  parameters  that  lead  to  the  lowest  AIC  value  should  be 
used.  The  statistical  parameters  that  maximize  the  log  likelihood  function  given  by  Eq.  (14) 
minimize  AIC. 


2.3.4. 1  Reference  S-N  curve  for  aluminum  24S-T3 


The  fatigue  data  of  93  unnotched  specimen  extracted  from  NACA-TN-2324,  Table  2  is  tabulated 
in  the  Excel  file:  naca-tn-24S-T3-(Raw  Data).xlsx,  worksheet  “Kt=l  data”.  Of  the  93  specimens, 
10  were  runouts  and  11  were  disqualified  tests1.  Fifty  three  (53)  data  points  were  used  for 
determining  the  S-N  curve  of  the  material  (29  specimens  tested  with  Smax  >  1.06  Sy,  where  Smax 
is  the  maximum  test  stress  in  the  load  cycle  and  Sy  is  the  average  tensile  yield  strength  of  the 
material  in  the  grain  direction  shown  in  Table  1,  were  excluded).  The  median  S-N  curve  of  the 
material  and  the  two  theoretical  design  curves  based  on  the  parameters  of  the  statistical  model  Si 
are  shown  in  Figure  3. 
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Reference  S-N  Curve  -  RFL  Normal  (53-data  points) 
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Figure  3:  Median  S-N  curve  for  aluminum  24S-T3  using  the  parameters  from  statistical  model  1. 


2.3.4.2  Reference  S-N  curve  for  aluminum  75S-T6 


1  Most  disqualified  tests  had  failures  that  occurred  outside  the  critical  section  (the  area  Vz-in  at  either  side  of  the  line 
of  minimum  cross-section).  A  smaller  number  of  failures  occurred  inside  the  grips  or  in  a  defect  (flaw  or  pit). 


Page  15 


ESRD,  Inc. 


SBIR  II  Option  3  -  Task  6 


N68335-1 6-C-0088 


The  fatigue  data  of  86  unnotched  specimen  extracted  from  NACA-TN-2324,  Table  3  is  tabulated 
in  the  Excel  file:  naca-tn-75S-T6-(Raw  Datal.xlsx,  worksheet  “Kt=l  data”.  Of  the  86  specimens, 
12  were  runouts  and  3  were  disqualified  tests.  Eighty  three  (83)  data  points  were  used  for 
determining  the  S-N  curve  of  the  material  (no  specimens  were  tested  with  Smax  >  1.06  Sy). 

The  median  S-N  curve  of  the  material  and  the  two  theoretical  design  curves  based  on  the 
parameters  of  the  statistical  model  Si  are  shown  in  Figure  4. 
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Figure  4:  Median  S-N  curve  for  aluminum  75S-T6  using  the  parameters  from  statistical  model  1. 


2.3.4.3  Reference  S-N  curve  for  SAE  4130  steel 


The  fatigue  data  of  59  unnotched  specimen  extracted  from  NACA-TN-2324,  Table  4  is  tabulated 
in  the  Excel  file:  naca-tn-SAE-4130-(Raw  Datal.xlsx,  worksheet  “Kt=l  data”.  Of  the  59 
specimens,  16  were  runouts  and  6  were  disqualified  tests.  Forty  Eight  (48)  data  points  were  used 
for  determining  the  S-N  curve  of  the  material  (5  specimens  tested  with  Smax  >  1.06  Sy,  where 
Smax  is  the  maximum  test  stress  in  the  load  cycle  and  Sy  is  the  average  tensile  yield  strength  of 
the  material  in  the  grain  direction  shown  in  Table  1,  were  excluded).  The  median  S-N  curve  of 
the  material  and  the  two  theoretical  design  curves  (P-quantiles)  based  on  the  parameters  of  the 
statistical  model  Si  are  shown  in  Figure  5. 
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Figure  5:  Median  S-N  curve  for  SAE-4130  steel  using  the  parameters  from  statistical  model  1. 

2.4  Generalization  for  Notches  -  Calibration  of  Predictors 


Predictors  of  damage  accumulation  are  phenomenological  models  constructed  for  the  purpose  of 
generalizing  sets  of  experimental  data  obtained  in  fatigue  tests  of  mechanical  or  structural 
components  subjected  to  cyclic  loading.  Predictors  are  scalars  that  depend  on  the  solution  of  a 
problem  of  elasticity,  or  more  generally,  a  problem  of  continuum  mechanics. 

The  predictors  used  in  classical  high  cycle  fatigue  life  estimation  depend  on  the  maximum 
normal  stress  amax  and  the  stress  amplitude  ca  defined  by 


cr  = 


^~max  *Tnin  


1  -R 


Eq.  23 


where  R,  called  the  cycle  ratio,  is  the  ratio  of  the  minimum  principal  stress  amjn  to  the  maximum 
principal  stress  at  the  location  of  the  maximum  principal  stress.  When  the  linear  theory  of 
elasticity  is  applicable  then  it  is  also  the  ratio  of  the  minimum  applied  load  to  the  maximum 
applied  load.  A  family  of  predictors  commonly  used  in  engineering  practice  for  uniaxial  stress  is: 

F  =  &eq  =  CTmax^a  ,  0  <  ^  <  1  Eq.  24 


Where  aeq  is  called  the  equivalent  stress.  Using  Eq.  (23): 

\-R 


F  ^"eq  ^"max 


0  <  q  <  1 


Eq.  25 


The  special  case  of  q  =  1/2  is  one  of  the  predictors  proposed  by  Smith,  Watson  and  Topper  in  [5] 
and  used  in  this  work  for  determining  the  S-N  curve  of  the  materials. 
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Predictors  are  formulated  with  the  objective  to  generalize  results  of  fatigue  tests  performed  on 
coupons  subjected  to  constant  or  smoothly  varying  uniaxial  or,  less  frequently,  biaxial  stress 
conditions  to  arbitrary  tri-axial  stress  conditions. 

Specifically,  we  investigated  the  performance  of  a  family  of  predictors  defined  as  follows: 


F  =  Ga=±-\(ccll+(\-a)cj)dV 

V  a, 


(l-R) 


1/2 


0<«<1 


Eq.  26 


where  I\  is  the  first  stress  invariant,  a  is  the  von  Mises  stress,  a  is  a  material-dependent 
parameter  determined  by  calibration  and  Vc  is  the  volume  of  the  domain  of  integration  defined  by 

nc  =  { x  I  cr,  >p(Tpeak  )  Eq.  27 


Where  3c  is  the  position  vector,  ai  is  the  first  principal  stress,  0  <  apcak,  is  the  maximum  stress  at 
the  stress  raiser,  J3  is  an  experimentally-determined  variable  that  depends  on  the  material,  the 
definition  of  Ga  and  the  geometric  features  of  the  stress  raiser. 

The  effect  of  the  geometric  features  of  the  stress  raisers  is  quantified  by  the  highly  stressed 
volume  HSV  which  is  defined  by: 

HSV  =  f  y(x)dV  where  y  =  f  WhCn  >  7  (7peak  Eq.  28 

j*  [0  otherwise 


The  domain  of  integration  Qp  is  the  neighborhood  of  a  stress  raiser.  Denoting  the  position  vector 
of  the  stress  raiser  by  x0  we  have  Q.p  =  {jc  l|jc-3c0|  <  p]  where  p  is  chosen  large  enough  to 

include  all  points  such  that  cr,  (x)  >  ycrpeak  and  small  enough  to  include  only  one  stress  raiser  or 
one  group  of  closely  spaced  stress  raisers. 

The  parameter  y  is  independent  of  material  properties.  Its  role  is  to  define  the  highly  stressed 
volume  HSV  which  depends  on  the  stress  field  in  the  neighborhood  of  stress  raisers.  In  a  separate 
investigation  it  was  established  that  the  predicted  number  of  fatigue  cycles  is  not  sensitive  to  the 
choice  of  y  and  y  =  0.85  is  a  reasonable  choice.  The  same  value  was  used  in  [1],  Henceforth  we 
will  also  use  y  =  0.85  in  this  report. 

An  example  of  the  highly  stressed  volume  for  a  circular  cutout  using  y  =  0.85  is  shown  in  Figure 

6. 
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Figure  6:  Example  of  the  stress  distribution  around  a  circular  cutout  in  a  test  specimen  (left)  and  the 
corresponding  size  of  the  HSV  computed  using  Eq.  (28)  with  y=  0.85  (right). 

In  the  special  case  of  uniaxial  stress  Ga  is  the  equivalent  uniaxial  stress.  The  integrand  is  the 
generalization  of  uniaxial  stress  to  triaxial  stress  conditions  considered  in  this  work.  Equations 
(26)  to  (28)  define  the  functional  form  of  the  predictor.  The  parameters  a  and  (i  are  material- 
dependent  and  thus  have  to  be  estimated  through  calibration. 

Remark:  To  the  authors’  knowledge,  using  a  linear  combination  of  the  first  stress  invariant  and 
the  octahedral  shear  stress  for  the  prediction  of  fatigue  life  under  tri-axial  loading  conditions  was 
first  proposed  by  Sines  in  [8].  There  are  important  differences  between  the  predictor  described 
by  Sines  and  the  predictor  proposed  herein:  (a)  Sines’  predictor  is  a  pointwise  quantity  whereas 
the  predictor  described  herein  is  an  integral  average  over  a  solution-dependent  domain,  (b)  Sines 
assumed  that  1\  is  determined  by  the  static  stress  and  the  octahedral  shear  stress  is  determined  by 
the  alternating  stress.  No  such  assumption  is  made  here:  Interaction  between  the  mean  and 
alternating  stresses  is  taken  into  account  through  the  equivalent  stress  aeq. 

2.4.1  Calibration  of  the  predictor 

The  calibration  of  predictors  is  based  on  the  test  records  shown  in  Table  2.  In  the  column  labeled 
Actual  Kt,  the  actual  stress  concentration  factors  determined  from  finite  element  analysis, 
defined  as  the  ratio  of  maximum  (peak)  stress  to  the  average  stress  in  the  test  section,  are  shown. 
In  the  column  labeled  HSA  the  highly  stressed  area  is  shown.  This  area,  multiplied  by  the 
thickness,  approximates  the  highly  stressed  volume  defined  in  Eq.  (28)  for  y  =  0.85.  In  the  last  3 
columns  the  specimen  count  for  each  material  is  listed  with  and  without  runouts.  The  number  of 
specimens  without  counting  runouts  is  in  brackets.  This  count  does  not  include  specimens  that 
were  tested  at  stress  levels  higher  than  1.06  Sy  and  specimens  that  were  disregarded  by  the 
investigators  because  of  flaws,  buckling,  or  failure  in  the  grips  as  explained  in  Section  2.3.4.  The 
test  records  included  in  the  count  are  called  qualified  test  records. 
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Each  specimen  type  is  identified  by  the  first  column  entry  in  Table  2.  For  each  specimen  there 
are  a  set  of  mk  qualified  test  records;  ( S(JJt ,  /V„  R,),  /  =  1,2,...,  mk  where  S<n  test  is  the  average 

stress  in  the  gauge  section  which  is  defined  as  the  minimum  width  times  the  thickness.  All  the 
specimen  types  and  the  corresponding  finite  element  mesh  for  each  are  shown  in  Appendix  3. 

The  highly  stressed  volume  HSV  depends  on  the  geometry  and  is  independent  of  the  magnitude 
of  loading.  Therefore  the  HSV  has  to  be  computed  for  each  specimen  type  only  once.  In  planar 
problems,  when  the  thickness  to  notch  radius  ratio  is  small,  the  highly  stressed  area  HSA 
multiplied  by  the  thickness  is  a  very  close  approximation  to  HSV.  Some  of  the  specimen  types 
considered  in  this  investigation  do  not  meet  this  criterion.  However  a  separate  investigation 
showed  that  the  difference  between  the  highly  stressed  volume  and  the  highly  stressed  area 
multiplied  by  the  thickness  is  negligibly  small. 

The  workflow  for  determining  the  value  of  fJ>  for  a  given  value  of  a  is  shown  schematically  in 
Figure  7.  In  the  calibration  process  we  fix  0  <  a  <  1  at  an  arbitrary  value  and  estimate  [i,  by 
iteration  as  follows. 

•  For  a  given  specimen  type  ©  we  find  by  experiment  the  number  of  cycles  A,  at  which  the 
ith  notched  specimen  fails  ©. 

•  If  the  assumption  that  GJfj)  is  a  proper  generalization  of  ceq  is  correct,  then  the 
probability  distribution  of  G(p  ( [3i )  is  the  same  as  that  of  aeq  and  its  expected  value  is 
/r'dogio  Ni )  (D  where  n  is  the  calibrated  mean  defined  in  Eq.  (3). 

•  We  find  the  value  of  Smax  from  Eq.  (1)  ©  and  the  corresponding  value  of  [j,  by  iteration 
©  such  that 

G«(A)-^1(log10Ai)  =  0  i  =  1,2,..., mk  Eq.  29 

•  This  results  in  a  sequence  of  ft  values  ©  that  depend  on  (a)  the  statistical  model,  (b)  the 
qualified  calibration  records,  (c)  the  parameter  a  and  (d)  the  stress  field  in  the 
neighborhood  of  the  stress  raiser  of  a  specimen  type,  characterized  by  the  highly  stressed 
area  HSA. 

For  example,  for  a  =  0.6  and  the  test  record  highlighted  in  Figure  7,  the  value  of  ceq  is  35.02  ksi 
and  for  R  =  -1  the  corresponding  Smax  has  the  same  value  (see  Eq.  (1)).  It  was  found  by  iteration 
that  for  P  =  0.664,  5max  =  34.99  ksi  which  differs  from  the  target  value  of  35.02  by  less  than 
0.1%.  This  example  is  shown  in  Figure  8  for  specimen  type  2  of  75S-T6,  which  also  shows  the 
variation  of  p  with  a  for  the  highlighted  test  case. 
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O s™  =  1  J (a. I,  +  (1  ■ - a)a)dV,  Q.c  =  { *  |  <r,  >  0, <rpeal } 

vc  a 


O 


fi(a) 


Figure  7:  Workflow  for  the  calibration  of  the  predictor. 
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0.2 

0  688 

0.4 
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0.6 
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1.0 
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R 

N 
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s 
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p 

S 
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(ksi) 

28,000 

-1.00 

14,000 

53.73 

53.73 

0.848 

53,74 

24,000 

-1.00 

42,000 

39.93 

39,93 

0.662 

39.92 

21,000 

-LOO 

86,000 

35.02 

35.02 
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34.99 

18,000 

-LOO 

412,400 

29.75 

29.75 

0.656 

29.77 

16,000 

-LOO 

1,028,000 

28.42 

28.42 

0.740 

28,41 

Figure  8:  Example  of  calculation  of  ft  for  a  =  0.6  for  several  test  results  of  specimen  type  2  (left)  and  the 
variation  of  ft  with  a  for  the  highlighted  test  result  (right). 

The  value  of  /?  for  all  eight  specimen  types,  all  test  cases,  and  for  all  three  materials  are  available 
for  seven  values  of  «in  the  Excel  files:  BetaSummarv24S-T3.xlsx,  BetaSummary75S-T6.xlsx, 
BetaSummaryS AE-4 1 3Q.xlsx  (see  Appendix  4  for  details). 

The  calibration  of  the  predictor  involves  the  following  steps: 

1.  Compute  fi  values  corresponding  to  a  sequence  0  <  cq,  <*2,  .  .  .  ,  an  <  1  for  each 
specimen  available  for  calibration.  Eight  test  records  identified  by  line  items  2  to  8  and 
10  in  Table  2  were  used  for  calibration  of  /?  as  a  function  of  the  HSA.  Two  types  of  fit 
were  used  for  this  relationship: 
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a.  Using  linear  regression  to  compute  an  estimate  of  the  mean  values  of  /?,  denoted 
by  P ,  for  each  a  in  the  form: 

/?(«,. ,  HSA)  =  JUp=  aj  +  bj  log10  HSA,  j  =  1, 2, ,  n.  Eq.  30 

This  fitting  assumes  that  for  a  fixed  HSA,  logio  /3  is  normally  distributed  with 
standard  deviation  sp.  The  three  parameters  of  the  statistical  model  (a,  b,  sp)  were 
estimated  by  the  maximum  likelihood  method. 

b.  Logarithm  fitting  of  the  arithmetic  average  of  [j,  denoted  by  [javg,  for  each  a  in 
the  form: 

Pavg  (aj  ’  HSA)  =  dj  +  bj  In  HSA,  j  =  1, 2, . . . ,  n.  Eq.  31 

The  parameters  of  the  logarithmic  fit  and  the  coefficient  of  determination  were 
obtained  using  the  Trend  line  Option  in  Excel. 

An  example  of  the  relation  [i  -  HSA  for  Aluminum  24S-T3  (a  =  0.2)  is  shown  in  Figure 
9.  The  value  of  the  arithmetic  average  of  ft  for  the  8  test  records  characterized  by  the  HSA 
is  shown  next  to  the  diamond  symbol  labeled  “beta_avg”  in  the  graph;  the  logarithmic  fit 
is  shown  with  a  red  solid  line  characterized  by  R2  =  0.933  and  labeled  “Log(beta_avg)”; 
while  the  mean  value  of  /?  by  linear  regression  is  indicated  with  the  green  dotted  curve 
labeled  as  “beta_mean  logN”. 


Figure  9:  Example  of  the  relation  between  ft  and  the  highly  stress  area  HSA  for  24S-T3  ( a  =  0.2). 
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The  parameters  of  the  two  types  of  fitting  indicated  above  are  presented  in  Figure  10  to 
Figure  12  for  aluminum  alloys  24S-T3  and  75S-T6  and  for  steel  SAE  4130,  respectively, 
for  seven  values  of  a.  The  data  is  also  available  in  the  Excel  file: 
FittedBetaSummary.xlsx  and  in  the  macro  Excel  files  used  for  the  computations: 
Prediction  (24S-T3  SWT)-RFL-Model-l-Alpha-xx.xlsm:  Prediction  (75S-T6  SWTVRFL- 
Model- 1  - Alpha-xx.xlsm:  Prediction-S AE-4 1 30-RFL-Model- 1  -Alpha-xx.xlsm:  with  xx  = 
00,  02,  04,  05,  06,  08,  10  (see  Appendix  4  for  details). 
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Figure  10:  Linear  regression  and  logarithmic  fitting  parameters  for  7  values  of  a  -  Aluminum  24S-T3 


j 

75S-T6 

logic  P  is  normally  distributed 

Logarithmic  Fit  | 

ai 

bj 

Std.  Dev 

ai 

bi 

R2 

1 

0.0 

0.00050 

0.05275 

0.06330 

0.94656 

0.03440 

0.920 

2 

0.2 

0.00001 

0.05675 

0.06706 

0.93919 

0.03598 

0.936 

3 

0.4 

-0.00112 

0.06101 

0.07089 

0.93029 

0.03755 

0.951 

4 

0.5 

-0.00241 

0.06307 

0.07267 

0.92463 

0.03825 

0.957 

5 

0.6 

-0.00398 

0.06522 

0.07462 

0.91822 

0.03892 

0.961 

6 

0.8 

-0.00838 

0.06942 

0.07815 

0.90360 

0.04012 

0.963 

7 

1.0 

-0.01525 

0.07331 

0.08140 

0.88531 

0.04099 

0.952 

Figure  11:  Linear  regression  and  logarithmic  fitting  parameters  for  7  values  of  a  -  Aluminum  75S-T6 
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Figure  12:  Linear  regression  and  logarithmic  fitting  parameters  for  7  values  of  a- Steel  SAE  4130 
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2.  Specimen  #9  (Figure  13)  is  suitable  for  the  selection  of  the  optimal  value  of  a  because 
the  highly  stressed  area  ( HSA  =5.908xl0'7  in2)  associated  with  the  small  notch  radius  (R 
=  0.0035  in)  is  outside  the  range  of  the  calibration  specimens  (see  Table  2).  In  the 

calibration  process  we  fixed  a  and  selected  such  that  lies  on  the  S-N  curve,  see 

Figure  7.  Consequently  any  value  of  /?  or  pAVG  in  the  range  of  calibration  is  an 
interpolation  based  on  the  assumption  that  the  functional  form  of  Ga  is  correct  for  any  a. 
The  following  procedure  was  used  to  select  the  optimal  value  for  a: 


Figure  13:  Specimen  #9  -  Notch  detail. 

a.  Using  Eq.  30,  compute  an  estimate  of  the  mean  value  of 
/3(aj,HSA),  j  =  1,2,  ...,n  for  the  test  records  available  for  the  calibration  of  the 

predictor. 

b.  Repeat  step  (a)  using  Eq.  31  to  compute  an  estimate  of  the  arithmetic  average  of 
j3AVG(a j, HSA),  j  =  1,2,  ...,n  for  the  test  records  available  for  the  calibration  of 

the  predictor. 

The  computed  values  of  p  (column  heading  (/?  )iogN)  and  pAVG  (column  heading 

(/?  )avg)  for  specimen  #9,  for  7  values  of  a,  and  for  all  three  materials  are  shown 
in  Table  3. 


Table  3:  Computed  values  of  /?  (Eq.  30)  and  Pavg  (Eq.  31)  for  specimen  #9 
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c.  Using  Eq.  (26)  compute  given  /3{aj,HSA)  and  j3AVG(aj,HSA), for  each  of 

the  test  records  corresponding  to  specimen  #9  and  for  each  material.  The  index  i 
refers  to  a  test  record.  The  results  are  included  in  the  Excel  files:  RankingData- 

24S-T3-Specimen-9.xlsx. _ RankingData-75S-T6-Specimen-9.xlsx,  and 

RankingData-SAE-4130-Specimen-9.xlsx.  (see  Appendix  4  for  details). 

d.  Using  Eq.  (14)  compute  the  log  likelihood  function  LLj  with  substituted  for 


cr[‘J .  The  value  of  a}  that  maximizes  the  log  likelihood  function  is  the  best  choice 


for  statistical  model  S\  and  the  available  calibration  data. 

In  order  to  assess  the  relative  merit  of  choosing  aj  over  cq,  given  the  statistical 
model  Si  and  the  calibration  data  D,  we  compute  the  Bayes  factors  characterized 
by  the  parameter  6X  (where  6-  [A  \  A  2  s  jut  s/},  see  Section  2.3.1): 


Pr(DIG  §,)  L,W  (JJ/S^ 

BF> W 


Eq.  32 


The  results  are  shown  in  Table  4.  The  largest  Bayes  factor  indicates 

the  optimal  choice  for  a,  which  are  highlighted  in  red  in  the  table. 


Table  4:  Bayes  factors  for  Specimen  #9 


Material 

yS-fit 

Bayers  Factors  for  a 

0.0 

0.2 

0.4 

0.5 

0.6 

0.8 

1.0 

24S-T3 

(^)logN 

1.00 

6.30 

39.38 

96.75 

230.30 

366.98 

16.22 

(/Oavg 

1.00 

4.50 

4.23 

0.97 

0.03 

<  10A(-8) 

<  10A(-27) 

75S-T6 

(^)logN 

1.00 

1.62 

15.45 

61.88 

420.8 

17899 

8.00E+05 

(/Oavg 

1.00 

8.55 

230.4 

1252 

4477 

2792 

0.01 

cac  /lion 

(/?)logN 

1.00 

0.021 

4.9E-04 

<  10A(-5) 

<10A(-10) 

<10A(-26) 

<10A(-62) 

5AE  4130 

(jff)AVG 

1.00 

7.11E-14 

The  workflow  for  the  selection  of  the  optimal  value  of  a  described  in  Step  2  is 


summarized  in  Figure  14. 
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p(a2)  =  -0.0271  l  +  0.061781og10(/ZSL4) 


Hi 

Failure 

ceq  (P)foroL 

0,0 

0.2 

0,4 

0.5 

0.6 

0.8 

1.0 

3000 

0 

72.49 

71.43 

70.12 

69.29 

68.19 

65.60 

62.18 

30000 

0 

57.99 

57.14 

56.10 

55.44 

54.55 

5248 

49.74 

104000 

0 

46.04 

46.15 

45.31 

44.77 

44.06 

4239 

40.18 

379000 

0 

35.69 

35.17 

34.52 

34.11 

33.57 

3230 

30.61 

760000 

0 

3123 

30.77 

30.21 

29.85 

29.38 

2826 

26.79 

1096000 

0 

27.68 

27.47 

26.97 

26.65 

26.23 

25.23 

23.92 

12-S7Q0O 

0 

25.09 

24.73 

24.27 

23.99 

23.61 

2271 

21.52  . 

11390000 

0 

25.09 

24.73 

24.27 

23.99 

23.61 

2271 

21.52 

1122GOOO 

1 

23.87 

23.52 

23.09 

2281 

22.45 

2160 

20.47 

11201000 

1 

22.30 

21.98 

21.58 

2132 

20.98 

20.13 

19.13 

13420000 

1 

22.30 

21.98 

21.58 

2132 

20.98 

20.13 

19.13 

10000 

0 

5833 

57.48 

56.43 

55.76 

54.33 

5279 

50.04 

10000 

0 

5111 

50.36 

49.44 

4885 

48.-03 

46.25 

43.34 

30500 

0 

43.71 

43.07 

42.28 

4178 

41.12 

39.55 

37.49 

11400-0 

0 

34.72 

34.21 

33.58 

33.19 

32.66 

3142 

29.78 

523000 

0 

24.94 

24.57 

24.12 

23.84 

23.46 

2257 

21.39 

1047500-0 

1 

16.73 

16.48 

16.18 

15.99 

15.74 

15.14 

14.35 

Bayes  factors: 

1.00 

6.30 

39.38 

96.75 

230  30 

366.98 

16.22 

—  \{a2Iy+{\-a2)cr)dV 

vc  nc 


Q.c={x\<jx>j3(a2)crpeak\ 


(T  S 

i — * 

1 

to 

eg  max 

l  2  J 

Figure  14:  Workflow  for  the  selection  of  the  optimal  value  of  a. 

Remark:  The  calibration  of  predictors  can  be  considered  reliable  only  if  the  numerical  accuracy 
of  the  computed  data  is  guaranteed.  The  values  of  the  predictor  reported  herein  were  computed 
with  StressCheck2.  The  relative  errors,  estimated  by  p-extension,  were  shown  to  be  under  1%. 

2.4.2  Summary  of  predictions 


In  this  section  we  present  a  summary  of  the  predictions  for  the  optimal  value  of  a  for  each 
material  and  for  the  two  options  used  for  the  calibration  of  /?  as  a  function  of  the  HSA.  For  each 
material,  the  summary  includes  the  calibration  data  for  the  value  of  a  that  maximizes  the  log 
likelihood  function  when  (a)  the  logio  P  is  assumed  to  be  normally  distributed,  Eq.  (30),  and  (b) 
fitting  the  arithmetic  mean  of  [j  with  a  logarithmic  fit,  Eq.  (31).  The  prediction  of  B  for  specimen 
#9  is  shown  in  graphical  form  for  both  //-fits,  together  with  the  representation  of  the  fatigue  data 
on  the  reference  S-N  curve  of  the  material.  Finally  the  prediction  of  fatigue  life  based  on  the 
mean  S-N  curve,  and  the  theoretical  design  curves  associated  with  p  =  1/1000  and  p  =  1/100000 
are  presented  in  tabular  form. 

2.4.2.1  Aluminum  24S-T3 


Calibration  data  24S-T3:  Two  optimal  values  of  a  are  presented  for  the  calibration  of  p  as  a 
function  of  the  highly  stressed  area  (HSA).  As  shown  in  Figure  15,  a=  0.8  maximizes  the  log 
likelihood  function  when  the  logio  P  is  assumed  to  be  normally  distributed  -  Eq.  (30),  while  a  = 


2  StressCheck  is  a  trademark  of  Engineering  Software  Research  and  Development,  Inc. 


Page  26 


ESRD,  Inc. 


SBIR  II  Option  3  -  Task  6 


N68335-1 6-C-0088 


0.2  maximizes  the  log  likelihood  function  when  fitting  the  arithmetic  mean  of  [j  with  a 
logarithmic  fit  -  Eq.  (31)  -  as  shown  in  Figure  16. 


# 

Specimen 

Pavg 

Pf.t  1 

Pl-OgN 

2 

Open  Hole  Kt=2.11 

0.7020 

0.6976 

0.6999 

8 

Edge  Notch  r=0.760  in 

0.6602 

0.6973 

0.6995 

3 

Edge  Notch  Kt=2.17 

0.6437 

0.6246 

0.6085 

4 

Fillet  Notch  Kt=2.19 

0.6241 

0.5668 

0.5446 

10 

Edge  Notch  r=0.071  in 

0.4945 

0.5262 

0.5039 

5 

Edge  Notch  Kt=4.43 

0.5241 

0.5033 

0.4822 

7 

Edge  Notch  Kt=5.83 

0.4207 

0.4598 

0.4436 

6 

Fillet  Notch  Kt=4.83 

0.4166 

0.4103 

0.4034 

Figure  15:  Calibration  data  for  a  =  0.8  maximizes  the  log  likelihood  function  when  the  login  P 's  assumed  to 
be  normally  distributed.  Tabular  data  with  the  arithmetic  average  of  ft  for  each  specimen  type  under  the 
header  jiw  g,  the  corresponding  fitted  value  using  the  logarithmic  fit  of  Eq.  (31)  under  /?l  n,  and  the  beta 
value  from  the  linear  regression  fit  of  Eq.  (30)  under  Aluminum  24S-T3. 


# 

Specimen 

Pavg 

Pfit 

Pl-OgN 

2 

Open  Hole  Kt=2.11 

0.7314 

0.7518 

I  0.7549 

8 

Edge  Notch  r=0.760  in 

0.7322 

0.7515 

1  0.7545 

3 

Edge  Notch  Kt=2.17 

0.7068 

0.6842 

|  0.6712  ! 

4 

Fillet  Notch  Kt=2.19 

0.6641 

0.6306 

|  0.6114 

10 

Edge  Notch  r=0.071  in 

0.5883 

0.5930 

1  0.5728 

5 

Edge  Notch  Kt=4.43 

0.6072 

0.5718 

1  0.5520 

7 

Edge  Notch  Kt=5.83 

0.5023 

0.5315 

1  0.5147 

6 

Fillet  Notch  Kt=4.83 

0.4676 

0.4856 

0.4752 

Figure  16:  Calibration  data  for  a  =  0.2  maximizes  the  log  likelihood  function  when  fitting  the  arithmetic 
mean  of  /?  with  a  logarithmic  fit.  Tabular  data  with  the  arithmetic  average  of  P  lor  each  specimen  type  under 
the  header  Pwt,,  the  corresponding  fitted  value  using  the  logarithmic  fit  of  Eq.  (31)  under  PyU,  and  the  beta 
value  from  the  linear  regression  fit  of  Eq.  (30)  under  /?LogN.  Aluminum  24S-T3. 

Note  that  both  curve  fits  are  shown  in  each  figure  to  illustrate  how  close  they  are  to  each  other, 

even  though  only  one  is  optimal  in  the  sense  of  Eq.  (32)  for  a  given  a.  The  box  with  solid  lines 

over  the  equation  indicates  the  expression  of  /?  for  which  the  optimal  value  of  a  maximizes  the 

log  likelihood  function,  while  the  box  with  dotted  lines  over  the  equation  indicates  the  fit  done 

for  the  other  expression  of  (3  for  the  same  value  of  a.  The  tabular  data  show  the  values  of  the 

arithmetic  average  of  [3  (heading  /?Avg)  for  all  8  calibration  specimens,  the  values  of  [3  obtained 

from  Eq.  (31)  (heading  (3  rt)  and  the  values  (3  obtained  from  Eq.  (30)  (heading  (3  L<)gN)-  Each  of 

the  last  two  columns  in  the  table  is  boxed  with  the  same  color  and  line  type  as  the  corresponding 

equations  in  the  graph. 
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Prediction  of  B  for  specimen  #9  24S-T3:  With  the  /?-fit  obtained  from  calibration,  the  value  of  /? 
for  the  HSA  of  specimen  #9  is  predicted  for  each  optimal  value  of  a  as  shown  in  Figure  17  and 
Figure  18  for  a  =  0.8  and  a  =  0.2,  respectively  (see  also  Table  3).  Additionally,  the  pair  (aeq  (/?), 
N)  for  each  of  the  17  tests  available  for  specimen  #9  are  shown  in  the  S-N  curve  of  the  material. 
Note  that  for  each  optimal  value  of  a,  one  set  of  test  points  is  more  evenly  distributed  around  the 
mean  S-N  curve  than  the  other.  In  Figure  17  the  open  circles  are  closer  to  the  mean  S-N  curve 
( /? -fit  from  Eq.  30)  and  in  Figure  18  the  solid  triangles  are  closer  to  the  mean  S-N  curve  ( /7-fit 
from  Eq.  31). 


oeq(/?)  fora  =0.8 


N  (cycles) 


P-HSA  fora  =0.8 


- *** - 1 

Predit 

values 

:ted 

H 

✓  ** 

of/? 

4 

/ 

/ 

0.4 

-  03 

0.2 

- w— 

— O' 

-  * 

♦ 

o 

▲ 

beta_avg 
beta  mean  logN 
beta_LogN 
beta_predict 

Log.  (beta avg) 

no 

1.0E-07  1.0E-06  1.0E-05  1.0E-04  1.0E-03  1.0E-02  1.0E-01  1.0E+00 

HSA 


Figure  17:  Predicted  value  of  j3  for  specimen  #9  (right)  and  plot  of  the  failed  specimens  on  the  reference  S-N 
curve  when  aeq  is  computed  from  the  predicted  value  of  (left).  Aluminum  24S-T3  for  a  -  0.8. 
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Figure  18:  Predicted  value  of  ji  for  specimen  #9  (right)  and  plot  of  the  failed  specimens  on  the  reference  S-N 
curve  when  aeq  is  computed  from  the  predicted  value  of  (left).  Aluminum  24S-T3  for  a  -  0.2. 
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Life  Prediction  24S-T3:  The  prediction  of  the  number  of  cycles  to  failure  for  the  17  tests 
available  for  specimen  #9  are  shown  in  Table  5  for  a  =  0.8  and  in  Table  6  for  a  =  0.2, 
respectively.  The  second  column  shows  the  value  of  aeq  (j3)  for  each  test,  the  next  three 
prediction  columns  show  the  mean  prediction  based  on  the  S-N  curve  of  the  material,  and  the 
predictions  based  on  the  theoretical  design  curves  with  P  =  1/1,000  and  with  P  =  1/100,000  using 
statistical  model  l(see  Figure  2).  The  next  column  shows  the  actual  fatigue  life  recorded  during 
the  test  (runouts  are  highlighted  in  red),  and  the  last  column  in  the  table  shows  a  “yes”  if  the 
predicted  number  of  cycles  to  failure  using  the  theoretical  design  curve  with  P  =  1/1,000  is  larger 
than  the  actual  fatigue  life. 


Table  5:  Prediction  of  the  number  of  cycles  to  failure  for  specimen  #9.  Material  24S-T3,  a  =  0.8,  //-fit  Eq.  (30) 


# 

Seq  (ksi) 

Prediction  of  N 

Actual 
Fatigue  Life 

Over 

Prediction 

P=l/1000 

o 

o 

o 

-H 

Mean 

Design  curve 
P=l/1000 

Design  Curve 
P=l/100000 

1 

65.60 

8,479 

3,285 

2,281 

8,000 

no 

o 

o 

o 

o 

o 

2 

52.48 

17,319 

6,438 

4,426 

30,000 

no 

o 

o 

o 

3 

42.39 

38,661 

13,055 

8,788 

104,000 

no 

i 

-H 

4 

32.30 

151,652 

37,146 

23,590 

379,000 

no 

I 

p 

CO 

05 

o 

5 

28.26 

411,846 

68,526 

41,219 

760,000 

no 

T 

o 

o 

6 

25.23 

1,565,370 

125,482 

70,186 

1,096,000 

no 

7 

22.71 

37,023,383 

245,855 

123,428 

1,287,000 

no 

1 

_ 

1 

8 

22.71 

37,023,383 

245,855 

123,428 

11,390,000 

no 

9 

21.60 

100,000,000 

358,139 

166,792 

11,226,000 

no 

1 

1 

10 

20.18 

100,000,000 

649,572 

262,075 

11,281,000 

no 

-s 

<M 

o 

11 

20.18 

100,000,000 

649,572 

262,075 

13,420,000 

no 

o 

O 

12 

52.79 

16,974 

6,320 

4,346 

10,000 

no 

-H 

13 

46.25 

27,342 

9,694 

6,593 

18,000 

no 

i— i 

LO 

14 

39.55 

52,139 

16,718 

11,136 

38,500 

no 

-► 

rH 

15 

31.42 

181,364 

41,868 

26,344 

114,000 

no 

( 

| 

16 

22.57 

60,273,340 

256,979 

127,952 

523,000 

no 

17 

15.14 

100,000,000 

100,000,000 

8,716,786 

10,475,000 

no 

l 

1 
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Table  6:  Prediction  of  the  number  of  cycles  to  failure  for  specimen  #9.  Material  24S-T3,  a  =  0.2,  fi-Tit  Eq.  (31) 


# 

Seq  (ksi) 

Prediction  of  N 

Actual 
Fatigue  Life 

- 1 

Over 

Prediction 
P= 1/1000 

Mean 

Design  curve 
P=l/1000 

Design  Curve 
P=l/100000 

1 

67.62 

7,746 

3,009 

2,091 

8,000 

no 

2 

54.10 

15,621 

5,855 

4,033 

30,000 

no 

3 

43.69 

34,149 

11,750 

7,940 

104,000 

no 

4 

33.29 

126,054 

32,679 

20,947 

379,000 

no 

5 

29.13 

316,949 

59,145 

36,102 

760,000 

no 

6 

26.01 

1,011,267 

105,667 

60,466 

1,096,000 

no 

7 

23.41 

8,605,444 

199,689 

104,010 

1,287,000 

no 

8 

23.41 

8,605,444 

199,689 

104,010 

11,390,000 

no 

9 

22.26 

372,052,934 

283,655 

138,581 

11,226,000 

no 

10 

20.81 

100,000,000 

489,913 

212,407 

11,281,000 

no 

11 

20.81 

100,000,000 

489,913 

212,407 

13,420,000 

no 

12 

54.42 

15,316 

5,749 

3,961 

10,000 

no 

13 

47.67 

24,394 

8,767 

5,980 

18,000 

no 

14 

40.77 

45,592 

14,980 

10,026 

38,500 

no 

15 

32.38 

149,057 

36,711 

23,334 

114,000 

no 

16 

23.26 

10,794,393 

208,149 

107,651 

523,000 

no 

17 

15.61 

100,000,000 

100,000,000 

4,482,878 

10,475,000 

no 

o 
o 

o 

-H 

t"- 

o 
o 

o 

u 


* 


o 

o 

Q 

o 

© 

H 

CO 

CTi 

O 

O 


n 


A3 

U 

u 

O 


LJ 


CN 

Q 

© 

o 

-H 

t-H 

to 


2A.2.2  Aluminum  75S-T6 

Calibration  data  75S-T6:  Two  optimal  values  of  a  are  presented  for  the  calibration  of  (3  as  a 
function  of  the  highly  stressed  area  ( HSA ).  As  shown  in  Figure  19,  a=  1.0  maximizes  the  log 
likelihood  function  when  the  logio  P  is  assumed  to  be  normally  distributed  -  Eq.  (30),  while  a  = 
0.6  maximizes  the  log  likelihood  function  when  fitting  the  arithmetic  mean  of  ft  with  a 
logarithmic  fit  -  Eq.  (31)  -  as  shown  in  Figure  20. 


# 

Specimen 

Pavg 

PFIT 

PLogN 

2 

Open  Hole  Kt=2.11 

0.7626 

0.7402 

0.7448 

8 

Edge  Notch  r=0.760  in 

0.6957 

I  0.7398 

0.7443 

3 

Edge  Notch  Kt=2. 17 

0.6657 

1  0.6622 

0.6478 

4 

Fillet  Notch  Kt=2.19 

0.6424 

0.6003 

0.5800 

10 

Edge  Notch  r=0.071  in 

0.5501 

1  0.5570 

0.5367 

5 

Edge  Notch  Kt=4.43 

0.5365 

0.5325 

0.5137 

7 

Edge  Notch  Kt=5.83 

0.4728 

0.4861 

0.4728 

6 

Fillet  Notch  Kt=4.83 

0.4256 

L  0.4331 

0.4301 

Figure  19:  Calibration  data  for  a  =  1.0  maximizes  the  log  likelihood  function  when  the  logio /?  is  assumed  to 
be  normally  distributed.  Tabular  data  with  the  arithmetic  average  of  ft  for  each  specimen  type  under  the 
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header  /?Avg>  the  corresponding  fitted  value  using  the  logarithmic  fit  of  Eq.  (31)  under  Pm,  and  the  beta 
value  from  the  linear  regression  fit  of  Eq.  (30)  under  P\AK%.  Aluminum  75S-T6. 


# 

Specimen 

Pavg 

Pfit 

Pl.ogN 

2 

Open  Hole  Kt=2.11 

0.7822 

0.7804 

|  0.7865  1 

8 

Edge  Notch  r=0.760  in 

0.7470 

0.7801 

|  0.7861  1 

3 

Edge  Notch  Kt=2.17 

0.7152 

0.7064 

|  0.6948  1 

4 

Fillet  Notch  Kt=2.19 

0.6742 

0.6476 

|  0.6296  1 

10 

Edge  Notch  r=0.071  in 

0.6199 

0.6065 

|  0.5877  | 

5 

Edge  Notch  Kt=4.43 

0.6010 

0.5832 

|  0.5652  | 

7 

Edge  Notch  Kt=5.83 

0.5337 

0.5391 

|  0.5250  | 

6 

Fillet  Notch  Kt=4.83 

0.4591 

1  0.4826  | 

Figure  20:  Calibration  data  for  a  =  0.6  maximizes  the  log  likelihood  function  when  fitting  the  arithmetic 
mean  of  p  with  a  logarithmic  fit.  Tabular  data  with  the  arithmetic  average  of  P  for  each  specimen  type  under 
the  header  /?avg>  the  corresponding  fitted  value  using  the  logarithmic  fit  of  Eq.  (31)  under  P\w,  and  the  beta 
value  from  the  linear  regression  fit  of  Eq.  (30)  under  yffLogN.  Aluminum  75S-T6. 

As  before,  both  curve  fits  are  shown  in  each  figure  to  illustrate  how  close  they  are  to  each  other, 

even  though  only  one  is  optimal  in  the  sense  of  Eq.  (32)  for  a  given  a.  The  box  with  solid  lines 

over  the  equation  indicates  the  expression  of  /?  for  which  the  optimal  value  of  a  maximizes  the 

log  likelihood  function,  while  the  box  with  dotted  lines  over  the  equation  indicates  the  fit  done 

for  the  other  expression  of  [3  for  the  same  value  of  a.  The  tabular  data  show  the  values  of  the 

arithmetic  average  of  / 3  (heading  (3 avg)  for  all  8  calibration  specimens,  the  values  of  f3  obtained 

from  Eq.  (31)  (heading  (3 fit)  and  the  values  (3  obtained  from  Eq.  (30)  (heading  [3 L<>gN)-  Each  of 

the  last  two  columns  in  the  table  is  boxed  with  the  same  color  and  line  type  as  the  corresponding 

equations. 

Prediction  of  /?  for  specimen  #9  75S-T6:  With  the  /?- fit  obtained  from  calibration,  the  value  of  [3 
for  the  HSA  of  specimen  #9  is  predicted  for  each  optimal  value  of  a  as  shown  in  Figure  21  and 
Figure  22  for  a  =  1.0  and  a  =  0.6,  respectively  (see  also  Table  3).  Additionally,  the  pair  (aeq  (/?), 
N)  for  each  of  the  17  tests  available  for  specimen  #9  are  shown  in  the  S-N  curve  of  the  material. 
Note  that  for  each  optimal  value  of  a ,  one  set  of  test  points  is  more  evenly  distributed  around  the 
mean  than  the  other.  In  Figure  21  the  open  circles  are  closer  to  the  mean  S-N  curve  (/? -fit  from 
Eq.  30)  and  in  Figure  22  the  solid  triangles  are  closer  to  the  mean  S-N  curve  (/?-fit  from  Eq.  31). 
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Figure  21:  Predicted  value  of  p  for  specimen  #9  (right)  and  plot  of  the  failed  specimens  on  the  reference  S-N 
curve  when  ct1(|  is  computed  from  the  predicted  value  of  /?  (left).  Aluminum  75S-T6  for  a  =  1.0. 


°eq  I/O  for  a=  0.6 


v\  \ 

o 

- Median  Fatigue  Curve 

—  —  Design  Curve  P=  1/1000 

- Design  Curve  P=  1/100000 

▲  Specimen  #9  (beta_avg) 

o  Specimen  #9  (beta LogN) 
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Figure  22:  Predicted  value  of  fi  for  specimen  #9  (right)  and  plot  of  the  failed  specimens  on  the  reference  S-N 
curve  when  a,(|  is  computed  from  the  predicted  value  of  P  (left).  Aluminum  75S-T6  for  a  =  0.6. 


Life  Prediction75S-T6:  The  prediction  of  the  number  of  cycles  to  failure  for  the  17  tests 
available  for  specimen  #9  are  shown  in  Table  7  for  a  =  1.0  and  in  Table  8  for  a  =  0.6, 
respectively.  The  second  column  shows  the  value  of  aeq  (/?)  for  each  test,  the  next  three 
prediction  columns  show  the  mean  prediction  based  on  the  S-N  curve  of  the  material,  and  the 
predictions  based  on  the  theoretical  design  curves  with  P  =  1/1,000  and  with  P  =  1/100,000  using 
statistical  model  l(see  Figure  2).  The  next  column  shows  the  actual  fatigue  life  recorded  during 
the  test  (runouts  are  highlighted  in  red),  and  the  last  column  in  the  table  shows  a  “yes”  if  the 
predicted  number  of  cycles  to  failure  using  the  theoretical  design  curve  with  P  =  1/1,000  is  larger 
than  the  actual  fatigue  life.  Note  that  only  for  test  #8  in  Table  8  the  prediction  of  the  number  of 
cycles  to  failure  (6,899,204)  is  larger  than  the  actual  fatigue  life  (5,045,000). 
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Table  7:  Prediction  of  the  number  of  cycles  to  failure  for  specimen  #9.  Material  75S-T6,  a  =  1.0,  /?-fit  Eq.  (30) 


# 

Seq  (ksi) 

Prediction  of  N 

Actual 
Fatigue  Life 

- ^ 

Over 

Prediction 

P=l/1000 

1.007  ±0.001 

0006 

Mean 

Design  curve 
P= 1/1000 

Design  Curve 
P=l/ 100000 

1 

2 

80.32 

4,892 

2,211 

1,463 

11,000 

no 

59.67 

10,319 

4,557 

2,999 

31,500 

no 

i 

\ 

O 

3 

50.49 

17,034 

7,170 

4,668 

93,000 

no 

2 

M 

CO 

o 

o 

o 

© 

4 

5 

39.01 

46,937 

16,115 

10,096 

217,000 

no 

32.13 

166,847 

34,644 

20,287 

353,000 

no 

6 

7 

8 

29.83 

394,392 

49,379 

27,644 

482,000 

no 

rG 

CJ 

w 

O 

22.95 

100,000,000 

355,396 

122,821 

3,375,000 

no 

p- 

1 

O 

o 

.  -H 
© 
t-H 

to 

1— 1 

r 

20.65 

100,000,000 

2,963,171 

338,536 

5,045,000 

no 

9 

10 

19.97 

100,000,000 

26,565,547 

538,306 

11,299,000 

no 

18.36 

100,000,000 

100,000,000 

4,348,178 

14,874,000 

no 

11 

12 

64.91 

8,227 

3,679 

2,428 

6,300 

no 

55.11 

12,966 

5,623 

3,685 

16,000 

no 

13 

44.97 

25,559 

10,104 

6,493 

26,000 

no 

14 

35.82 

74,573 

22,048 

13,494 

39,500 

no 

15 

30.69 

269,378 

42,863 

24,460 

149,000 

no 

16 

17 

25.66 

100,000,000 

124,145 

58,649 

2,441,000 

no 

19.03 

100,000,000 

100,000,000 

1,347,182 

14,390,000 

no 

Table  8:  Prediction  of  the  number  of  cycles  to  failure  for  specimen  #9.  Material  75S-T6,  a  =  0.6,  /?- fit  Eq.  (31) 


# 

Seq  (ksi) 

Prediction  of  N 

Actual 
Fatigue  Life 

- 1 

Over 

Prediction 

P=l/1000 

Mean 

Design  curve 
P= 1/1000 

Design  Curve 
P=l/100000 

l 

78.88 

5,100 

2,305 

1,525 

11,000 

no 

2 

58.60 

10,854 

4,777 

3,141 

31,500 

no 

3 

49.58 

18,080 

7,551 

4,907 

93,000 

no 

4 

38.31 

51,363 

17,180 

10,717 

217,000 

no 

5 

31.55 

198,674 

37,616 

21,817 

353,000 

no 

6 

29.30 

529,988 

54,247 

29,958 

482,000 

no 

7 

22.54 

100,000,000 

449,293 

141,899 

3,375,000 

no 

8 

20.28 

100,000,000 

6,899,204 

427,885 

5,045,000 

yes 

9 

19.61 

100,000,000 

100,000,000 

725,382 

11,299,000 

no 

10 

18.03 

100,000,000 

100,000,000 

12,735,275 

14,874,000 

no 

11 

63.75 

8,627 

3,850 

2,539 

6,300 

no 

12 

54.12 

13,689 

5,906 

3,866 

16,000 

no 

13 

44.17 

27,395 

10,688 

6,851 

26,000 

no 

14 

35.18 

83,476 

23,658 

14,393 

39,500 

no 

15 

30.14 

339,824 

46,855 

26,423 

149,000 

no 

16 

25.20 

100,000,000 

142,696 

65,196 

2,441,000 

no 

17 

18.69 

100,000,000 

100,000,000 

2,221,527 

14,390,000 

no 
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2.4.23  Steel  SAE  4130 


Calibration  data  SAE  4130:  For  this  material,  one  value  of  a  was  found  to  be  optimal  for  the 
calibration  of  f3  as  a  function  of  the  highly  stressed  area  ( HSA ).  As  shown  in  Figure  23,  a=  0.0 
maximizes  the  log  likelihood  function  when  the  logio  (3  is  assumed  to  be  normally  distributed  - 
Eq.  (30)  -  and  when  fitting  the  arithmetic  mean  of  /3  with  a  logarithmic  fit  -  Eq.  (31). 


P-HSAfor  a  =0.0 


PAVG  =0.926+  0.04304  In  (HSA) 
R2  =0.892 

0.43 

0.4871 
£3 _ i _ 

B 

I 

Up  =  0.01361  +  0.07681og10(i/&4) 
5  =  0.0712 


♦  beta_avg 

—  —  betamean  logN 

Log.  (beta_avg) 

0.001 

HSA 


Specimen 

Pmg 

p™ 

PlogU 

Open  Hole  Kt=2.11 

0.7187 

0.7735 

0.7862 

Edge  Notch  r=0.760  in 

0.7795 

0.7732 

0.7857 

Edge  Notch  Kt=2. 17 

0.7041 

0.6917 

0.6793 

Fillet  Notch  Kt=2. 19 

0.6804 

0.6267 

0.6050 

Edge  Notch  r=0.071  in 

0.6424 

0.5812 

0.5578 

Edge  Notch  Kt=4.43 

0.5204 

0.5555 

0.5328 

Edge  Notch  Kt=5.83 

0.4878 

0.5067 

0.4884 

Fillet  Notch  Kt=4.83 

0.4263 

0.4511 

0.4423 

Figure  23:  Calibration  data  for  a  =  0.0  maximizes  the  log  likelihood  function  for  both  forms  of  /}  fitting. 
Tabular  data  with  the  arithmetic  average  of  ft  for  each  specimen  type  under  the  header  fiwt,,  the 
corresponding  fitted  value  using  the  logarithmic  fit  of  Eq.  (31)  under  /?Frr,  and  the  beta  value  from  the  linear 

regression  fit  of  Eq.  (30)  under  Steel  SAE  4130. 


Note  that  both  curve  fits  shown  in  the  figure  are  boxed  with  solid  lines  over  the  equation  since 
both  expressions  of  (3  are  optimal  for  the  same  value  of  a.  The  tabular  data  show  the  values  of 
the  arithmetic  average  of  /3  (heading  [3 avg)  for  all  8  calibration  specimens,  the  values  of  [3 
obtained  from  Eq.  (31)  (heading  /?fit)  and  the  values  f3  obtained  from  Eq.  (30)  (heading  /?u>gN)- 
Each  of  the  last  two  columns  in  the  table  is  boxed  with  the  same  color  and  line  type  as  the 
corresponding  equations. 

Prediction  of  [3  for  specimen  #9  SAE  4130:  With  the  /?-fit  obtained  from  calibration,  the  value  of 
f. 3  for  the  HSA  of  specimen  #9  is  predicted  for  the  optimal  value  of  a  as  shown  in  Figure  21  for  a 
=  0.0  (see  also  Table  3).  Additionally,  the  pair  (aeq  (/?),  N)  for  each  of  the  16  tests  available  for 
specimen  #9  are  shown  in  the  S-N  curve  of  the  material.  The  open  circles  were  determined  for 
the  f3  -fit  given  by  Eq.  30  and  the  solid  triangles  from  the  (3 -fit  from  Eq.  31. 
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P-HSAfor  a=  0.0 

Predicted 

^  * 

values  of  B 

Qrs 

/ 

1 

i 

05000 

- 

♦  beta_avg 

O  betaLogN 
—  —  beta  mean  LogN 

▲  beta_predict 
- Loo.  (beta  avo) 

0.0000001  0.000001  0.00001  0.0001  0.001  0.01  0.1  1 

HSA 

Figure  24:  Predicted  value  of  P  for  specimen  #9  (right)  and  plot  of  the  failed  specimens  on  the  reference  S-N 
curve  when  (jeq  is  computed  from  the  predicted  value  of  P  (left).  Steel  SAE  4130  for  a  =  0.0. 


Life  Prediction  SAE  4130:  The  prediction  of  the  number  of  cycles  to  failure  for  the  16  tests 
available  for  specimen  #9  are  shown  in  Table  9  when  //-fit  given  by  Eq.  30  and  in  Table  10 
when  the  //-fit  given  by  Eq.  31,  both  for  a  =  0.0.  The  second  column  shows  the  value  of  aeq  (//) 
for  each  test,  the  next  three  prediction  columns  show  the  mean  prediction  based  on  the  S-N  curve 
of  the  material,  and  the  predictions  based  on  the  theoretical  design  curves  with  P  =  1/1,000  and 
with  P  =  1/100,000  using  statistical  model  l(see  Figure  2).  The  next  column  shows  the  actual 
fatigue  life  recorded  during  the  test  (runouts  are  highlighted  in  red),  and  the  last  column  in  the 
table  shows  a  “yes”  if  the  predicted  number  of  cycles  to  failure  using  the  theoretical  design  curve 
with  P  =  1/1,000  is  larger  than  the  actual  fatigue  life.  Note  that  only  for  test  #6  in  Table  9  the 
prediction  of  the  number  of  cycles  to  failure  (330,840,339)  is  larger  than  the  actual  fatigue  life 
(3,193,000),  while  for  tests  #12  and  #14  in  Table  10  the  predicted  numbers  of  cycles  to  failure 
are  larger  than  the  actual  fatigue  life  (206.229  >  195,000  for  #12  and  824,410  >  389,000  for 
#14). 
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Table  9:  Prediction  of  the  number  of  cycles  to  failure  for  specimen  #9.  Steel  SAE  4130,  a  =  0.0,  ft- fit  Eq.  (30) 

# 

Seq  (ksi) 

Prediction  of  N 

Over 

Prediction 

P=l/1000 

007  ±  0.001 

)006 

Mean 

Design  curve 
P=l/1000 

Design  Curve 
P=l/100000 

Actual  Fatigue 

Life 

1 

79.71 

55,355 

7,763 

3,674 

124,000 

no 

© 

!L 

i _ > 

o 

2 

72.88 

95,664 

12,115 

5,631 

215,000 

no 

-H 

CO 

3 

63.77 

308,205 

26,130 

11,515 

379,000 

no 

t 

o 

o 

4 

54.66 

45,449,762 

83,187 

31,589 

919,000 

no 

o 

© 

5 

45.55 

100,000,000 

935,777 

178,161 

1,546,000 

no 

6 

41.00 

100,000,000 

330,840,339 

1,015,979 

3,193,000 

yes 

n 

— 1 

7 

36.44 

100,000,000 

100,000,000 

100,000,000 

14,394,000 

no 

8 

98.54 

20,948 

3,122 

1,498 

46,000 

no 

( 

/ 

o 

o 

9 

89.26 

31,662 

4,679 

2,239 

71,000 

no 

-s 

10 

70.57 

120,624 

14,395 

6,629 

182,000 

no 

o 

% 

o 

^•2 

11 

61.11 

542,964 

34,660 

14,853 

347,000 

no 

12 

56.34 

4,033,665 

64,014 

25,352 

195,000 

no 

r-H 

kO 

13 

56.34 

4,033,665 

64,014 

25,352 

874,000 

no 

-► 

14 

51.53 

100,000,000 

149,014 

50,494 

389,000 

no 

/ 

\ 

15 

51.53 

100,000,000 

149,014 

50,494 

12,508,000 

no 

16 

46.68 

100,000,000 

582,073 

133,160 

10,950,000 

no 

L 

l 

Table  10:  Prediction  of  the  number  of  cycles  to  failure  for  specimen  #9.  Steel  SAE  4130,  a  =  0.0 

\,P- fit  Eq.  (31 

) 

# 

Prediction  of  N 

Actual 
Fatigue  Life 

- 1 

Over 

Prediction 

P=l/1000 

.007  ±0.001 

1006 

Seq  (ksi) 

Mean 

Design  curve 
P=l/1000 

Design  Curve 
P=l/100000 

1 

70.89 

116,590 

14,044 

6,477 

124,000 

no 

o 

1 

o 

2 

64.81 

257,408 

23,600 

10,493 

215,000 

no 

iL 

\ 

-H 

CO 

05 

o 

Q 

3 

56.71 

3,090,802 

60,683 

24,224 

379,000 

no 

4 

48.61 

100,000,000 

306,757 

86,437 

919,000 

no 

© 

5 

40.51 

100,000,000 

100,000,000 

1,353,097 

1,546,000 

no 

j— 

"I 

6 

36.46 

100,000,000 

100,000,000 

100,000,000 

3,193,000 

no 

7 

32.41 

100,000,000 

100,000,000 

100,000,000 

14,394,000 

no 

8 

87.63 

34,427 

5,064 

2,420 

46,000 

no 

\ 

■8 

1  8 

o 

9 

79.38 

56,659 

7,920 

3,746 

71,000 

no 

10 

62.75 

374,889 

28,997 

12,654 

182,000 

no 

% 

o 

-H 

CO 

i-H 

to 

11 

54.35 

186,773,493 

87,668 

32,989 

347,000 

no 

12 

50.10 

100,000,000 

206,229 

64,720 

195,000 

yes 

13 

50.10 

100,000,000 

206,229 

64,720 

874,000 

no 

-► 

♦ 

14 

45.83 

100,000,000 

824,410 

165,254 

389,000 

yes 

L 

J 

15 

45.83 

100,000,000 

824,410 

165,254 

12,508,000 

no 

16 

41.51 

100,000,000 

47,955,657 

776,663 

10,950,000 

no 

l 

J 
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2.4.3  Validation  of  the  predictions 

The  statement:  “this  model  has  been  validated”  should  be  understood  to  mean  that  one  has  a  high 
degree  of  belief,  supported  by  experimental  evidence,  that  the  model  in  question  will  perform  its 
predictive  function  reliably.  This  statement  is  incomplete  however  without  a  qualifying  statement 
concerning  the  range  of  parameters  and  range  of  data  for  which  one  believes  the  model  has  been 
validated. 

We  considered  the  formulation  and  calibration  of  a  model  characterized  by  two  parameters  a  and 
(3.  We  used  a  subset  of  the  available  data,  the  training  set,  to  calibrate  the  parameter  /?  for  fixed 
values  of  a.  We  then  used  the  remaining  data,  the  validation  set,  to  determine  optimal  values  for 
a.  In  the  calibration  process  the  model  is  assumed  to  be  correct  and  /3i,  corresponding  to  the 
training  data,  are  determined  by  projecting  each  data  point  onto  the  S-N  curve  by  means  of  Eq. 
(29).  We  calibrated  the  model  for  each  a.  We  then  considered  ranking  the  models  with  respect  to 
the  validation  set  based  on  specimen  #9. 

As  shown  in  Figure  17  for  aluminum  24S-T3  (a  =  0.8),  the  predicted  values  of  (3  fall  outside  the 
range  of  the  training  set  (the  same  is  true  for  the  other  materials  and  for  different  values  of  a). 
We  then  updated  the  training  set  to  include  the  validation  set  and  obtained  a  new  calibration 
curve  p  -  HSA  for  each  material  and  for  the  optimal  values  of  a  corresponding  to  the  /Tfits  using 
Eq.  (30)  and  Eq.  (31).  The  updated  calibration  curves  are  shown  in  Figure  25  to  Figure  29. 
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Figure  25:  Updated  /?  calibration  for  a  =0.2  including  all  9  specimen  types.  Aluminum  24S-T3. 
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Figure  26:  Updated  P  calibration  for  a  -  0.8  including  all  9  specimen  types.  Aluminum  24S-T3. 
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Figure  27:  Updated  P  calibration  for  a  -  0.6  including  all  9  specimen  types.  Aluminum  75S-T6. 
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Figure  28:  Updated  P  calibration  for  a  -  1.0  including  all  9  specimen  types.  Aluminum  75S-T6. 
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Figure  29:  Updated  /? calibration  for  a=  0.0  including  all  9  specimen  types.  Steel  SAE  4130. 

The  data  is  also  available  in  the  Excel  files:  Calibration- All- ( 24S -T3  -S  WT)-RFL-Model- 1  - 
Alpha-02.xlsm;  Calibration- A11-124S-T3-S  WT)-RFL-Model- 1  -Alpha-08.xlsm;  Calibration- All- 
(75S-T6-SWT)-RFL-Model-l-Alpha-06.xlsm:  Calibration- All-(75S-T6-SWT)-RFL-Model- 1- 
Alpha- 1 0.xlsm;  Calibration- All-SAE-4 1 30-RFL- Model- 1  - Alpha-OO.xlsm  (see  Appendix  4  for 
details). 

We  now  state  that 

1.  The  model  has  been  validated.  Using  statistical  model  1  for  characterizing  the  S-N  data 
and  two  models  to  characterize  the  -  HSA  relation  for  each  material,  we  found  the 
optimal  values  of  a  that  maximizes  the  Bayes  factors  as  shown  in  Table  11  (extracted 
from  Table  4). 

Table  11:  Optimal  value  of  a  for  each  material  and  /?-fit 


Material 

Optimal  Value  of  a 

(ft  )logN 

(ft)  AVG 

24S-T3 

0.80 

0.20 

75S-T6 

1.00 

0.60 

SAE  4130 

0.00 

0.00 

2.  The  ranges  of  parameters  and  data  for  which  validation  is  claimed  are  approximately 
ICC7  <  HSA  <  10” 1  (in2),  104  <  N  <  107,  p  -  2s  <  loglO  N  <  p  +  2s  where  s  is  the  standard 
deviation. 

2.4.4  Formulation  of  design  rules 

Design  rules  for  structural  and  mechanical  components  subjected  cyclic  loading  are  stated  in  the 
form  of  a  predictor  and  its  allowable  value. 
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The  predictor  considered  herein  depends  on  the  material  properties  (because  the  parameter  a 
depends  on  the  material),  its  allowable  value  depends  on  both  the  material  properties  and  the 
number  of  load  cycles  that  the  component  is  expected  to  endure.  Incorporated  in  the  allowable 
value  is  a  factor  of  safety  chosen  to  account  for  both  aleatory  and  epistemic  uncertainties.  We 
distinguish  between  two  kinds  of  epistemic  uncertainties: 

1.  Model  form  uncertainties.  These  are  uncertainties  associated  with  the  assumptions 
incorporated  in  the  mathematical  problem,  the  definition  of  the  predictor  and  the 
definition  of  the  statistical  model. 

2.  Errors  in  the  numerical  solution.  It  must  be  verified  that  errors  in  the  numerical  solution 
of  the  mathematical  model  are  within  tight  tolerances  at  all  times.  Otherwise  model  form 
uncertainties  and  numerical  errors  cannot  be  separated,  making  proper  assessment  and 
ranking  of  mathematical  models  impossible.  Solution  verification  is  an  essential  technical 
requirement  of  numerical  simulation  [9]. 

We  have  considered  constant  amplitude  loading  only.  Additional  epistemic  uncertainties  of  the 
first  kind  are  associated  with  the  choice  of  cumulative  fatigue  damage  models  in  variable 
amplitude  loading.  Ideally,  epistemic  uncertainties  would  be  negligibly  small  in  comparison  with 
aleatory  uncertainties. 

2.5  Uncertainty  Quantification 

The  aleatory  uncertainties  in  (i  are  correlated  with  the  uncertainties  in  the  S-N  data.  Since  the  S- 
N  curve  represents  the  mean  fatigue  life  of  the  material,  the  mean  value  of  /?  for  a  given  Ga  is 
one  possible  generalization  for  the  notched  specimens.  We  utilize  the  generalization  of  the  mean 
value  of  [3  given  by  Eq.  (30)  as  described  in  Section  2.4.1,  which  assumes  that  for  a  fixed  HSA, 
logio  J3  is  normally  distributed  with  standard  deviation  sp.  Next  we  consider  the  mean  value  of  J3 
plus  one  standard  deviation  to  quantify  the  uncertainty  in  the  predicted  fatigue  life  of  Specimen 
#9. 

The  model-form  uncertainty  is  associated  with  the  selection  of  Ga ■  The  selection  of  Ga  depends 
on  the  model  used  for  fitting  the  available  data  during  calibration  for  determining  the  relation  fi  - 
HSA.  We  utilize  the  optimal  values  of  a  obtained  for  the  two  generalizations  of  the  mean  value 
of  P  given  by  Eq.  (30)  and  (31)  to  compute  two  values  of  Ga  for  the  prediction  of  fatigue  life  of 
Specimen  #  9. 

Aluminum  24S-T3.  Consider  the  quantification  of  the  aleatory  and  model-form  uncertainties 
associated  with  the  prediction  of  fatigue  life  of  specimen  #9. 
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The  predicted  value  of  /?  for  each  optimal  value  of  a  identified  in  Section  2.4.1  and  2.4.2  was 
used  to  quantify  one  aspect  of  the  model-form  uncertainty  in  the  prediction  of  the  mean  number 
of  cycles  to  failure  of  test  case  4  for  which  the  actual  fatigue  life  during  the  test  was  379,000 
cycles.  As  shown  in  Table  5  and  Table  6,  the  predicted  mean  number  of  cycles  to  failure  for  the 
model  based  on  Go.s  (corresponding  to  the  /?- fit  given  by  Eq.  30)  is  151,652,  while  for  the  model 
based  on  G0.2  (corresponding  to  the  /?-fit  given  by  Eq.  31)  is  126,054. 

Quantification  of  uncertainty  associated  with  the  choice  of  fitting  algorithm  for  / 3  is  summarized 
in  Figure  30  for  test  case  #4.  For  all  other  test  cases  the  effect  of  model-form  uncertainty  in  the 
prediction  of  the  number  of  cycles  to  failure  is  realized  by  comparing  the  ‘Prediction  of  N’ 
columns  of  Table  5  and  Table  6. 


Specimen  # 9  -  Test  #4 
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Figure  30:  Quantification  of  uncertainty  for  associated  with  the  selection  of  P  for  Specimen  #9  -  Test  #4. 

Aluminum  24S-T3. 


The  quantification  of  aleatory  uncertainty  is  illustrated  with  the  determination  of  ft  for  specimen 
#9  considering  the  mean  of  the  /?-fit  given  by  Eq.  (30)  and  the  mean  plus  one  standard  deviation. 
As  shown  in  Figure  31,  the  predicted  values  of  ft  for  each  case  are  0.3163  and  0.3925  for  the 
mean  curve  fit  and  the  mean  plus  one  standard  deviation,  respectively. 
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Figure  31:  Predicted  values  of  /}  for  specimen  #9  for  the  model  characterized  by  a  =  0.8  from  the  mean  /?- 
HAS  curve  fit  and  the  mean  plus  one  standard  deviation.  Aluminum  24S-T3. 


With  the  predicted  values  of  /?  we  now  compute  the  Ga  for  test  case  4  and  the  corresponding 
estimated  number  the  cycles  to  failure  using  the  mean  S-N  curve  of  the  material: 


lr  ( \  —  R  V/2 

%  =  Go.s  =  —  j  (0‘8 /1  +  °-2  ^ )dV  —Z~  with  CLc  =  {x\<71>/3  CTpeak } 

vcnr  V  ^  J 


For  p  =  0.316  — >  (Teq  =  32.30  ksi 
For  fi  =  0.392  ->  cre?  =  41.47  ksi 

Figure  32  shows  the  effect  of  the  aleatory  uncertainty  in  the  prediction  of  fatigue  life  assuming 
that  the  model  is  correct.  By  model  we  understand  the  entire  process  by  which  predictions  are 
made. 
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Figure  32:  Quantification  of  aleatory  uncertainty  for  Specimen  #9  -  Test  #4.  Aluminum  24S-T3. 


Steel  SAE  4130.  Consider  the  quantification  of  the  model-form  uncertainty  associated  with  the 
prediction  of  the  mean  number  of  cycles  to  failure  for  specimen  #9  of  steel.  In  this  case  we  look 
at  the  prediction  for  test  cases  #2  and  #10  shown  in  Table  9  and  Table  10  when,  for  the  same 
optimal  values  of  a=  0,  the  value  of  (i  is  given  by  Eq.  (30)  or  (31)  in  the  computation  of  Ga= o. 


These  specimens  were  tested  a  different  stress  levels  (Stest  =  40  ksi  and  Stest  =  50  ksi  for  #2  and 
#10,  respectively)  and  load  ratios  (R  =  — 1.0  and  R  =  —0.2,  for  #2  and  #10  respectively),  but  they 
have  very  similar  equivalent  stress,  and  therefore  the  expected  life  should  also  be  similar.  The 
effect  of  the  model-form  uncertainty  in  the  prediction  is  summarized  in  Figure  33. 
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Figure  33:  Quantification  of  model-form  uncertainty  for  Specimen  #9  -  Tests  #2  and  #10.  Steel  SAE  4130. 

2.5.1  Frequency  plots 

In  Section  2.4.1,  it  was  assumed  that  for  a  fixed  value  of  the  highly  stressed  area  ( HSA ),  logio  /?  is 
normally  distributed  and  that  the  mean  (jj.fi)  is  a  linear  function  of  logio  USA  with  constant 
standard  deviation  sp.  In  order  to  check  this  assumption,  a  frequency  plot  was  constructed  for  the 
optimal  value  of  a  for  each  material. 

Figure  42  shows  the  calibration  curve  of  /?  for  aluminum  24S-T3  for  a  =  0.8.  The  symbols 
correspond  to  the  computed  values  of  /?  for  each  test  case  of  the  8  specimen  types  used  for 
calibration  shown  in  Table  2.  The  dotted  line  is  the  estimate  of  the  mean  value  of  /?  as  a  function 
of  the  HSA  determined  from  Eq.  (30): 

Up  =  -0.04 184  +  0.07355  log10(//&4),  with  sp  =  0.094 

The  Z-score  measures  how  many  standard  deviations  each  point  f ?,  is  from  the  mean  /up.  The  Z- 
score  is  defined  as: 


Z.  = 


logi  oPi-Mp 


Eq.  33 
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The  groupings  of  the  Z-scores  in  bins  form  the  frequency  plot.  The  data  is  available  in  the  Excel 
file:  BetaZ-Score-Alpha-0.8(24S-T3).xlsx  and  summarized  in  Figure  35  to  Figure  37  for  each 
specimen  type  and  for  all  specimens  combined.  Since  specimens  2  and  8  have  practically  the 
same  HSA  (see  Table  2),  a  single  frequency  plot  was  prepared  for  both.  The  results  clearly  show 
that  the  log-normal  distribution  assumption  is  supported  by  the  resulting  distribution  in  the 
frequency  plots. 


Figure  34:  Calibration  curve  for  /},  aluminum  24S-T3  and  a  =  0.8.  Symbols  are  the  computed  values  of  P  for 
each  specimen  type  (Section  2.4.1),  and  the  dotted  line  is  the  fit  of  the  mean  value  of  /?. 


Figure  35:  Frequency  plot  for  specimen  types  2  &  8  combined,  3  and  4.  Aluminum  24S-T3. 
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Figure  36:  Frequency  plot  for  specimen  types  5, 6  and  7.  Aluminum  24S-T3. 


Figure  37:  Frequency  plot  for  specimen  type  10  (left)  and  for  all  the  specimens  combined  (right).  Aluminum 

24S-T3. 

The  combined  frequency  plots  for  all  the  specimen  types  of  aluminum  75S-T6  ( a=  1.0)  and  steel 
SAE-4130  (a=  0.0)  are  shown  in  Figure  38.  The  complete  data  set  for  each  material  is  available 
in  the  Excel  files:  BetaZ-Score- Alpha-1. 0(75S-T6).xlsx  and  BetaZ-Score-Alpha-0.0(SAE- 
4130).xlsx.  respectively. 
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Figure  38:  Frequency  plots  for  all  the  specimens  combined.  Aluminum  75S-T6  and  Steel  SAE-4130. 


Finally,  a  quantile  of  input  sample  v.  quantile  of  lognormal  distribution  for  all  9  specimen  types 
of  aluminum  24S-T3  for  a  =  0.8  (a  QQ  plot)  is  shown  in  Figure  39. 


QQ  Plot  of  normalized  logic  P  versus  Standard  Normal 


Standard  Normal  Quantiles 

Figure  39:  QQ  plot  for  all  9  specimen  types  of  aluminum  24S-T3  for  a  -  0.8. 
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3  Summary  and  Conclusions 

The  focus  of  activities  in  the  Option  3  Period  of  this  contract  was  to  expand  the  work  initiated 
during  the  Option  1  &  2  periods  [10,  11]  for  the  support  of  flaw -tolerant  design,  certification  and 
sustainment  based  on  a  reliable  predictive  capability  for  the  estimation  of  the  expected  life  of 
metallic  airframe  components,  given  part  specifications,  material  properties,  initial  flaws  and 
design  load  spectra. 

During  the  Option  2  Period  the  focus  was  on  providing  a  framework  and  recommendations  for 
the  development  of  stress-based  damage  accumulation  models  including  the  extension  of  the 
models  for  multiaxial  loading.  Stress-based  damage  accumulation  models  using  the  integral 
average  of  a  function  of  stress  over  a  solution-dependent  volume  characterized  by  a  single 
parameter  {(3)  provided  a  basis  for  ranking  in-service  damage  from  the  point  of  view  of  fatigue 
criticality. 

During  the  Option  2  Period  we  found  that  the  best  results  were  obtained  when  the  highly  stressed 
volume  ( HSV)  of  material  around  a  flaw  was  correlated  with  /?,  allowing  the  generalization  of  the 
calibration  work  performed  during  Option  1  Period  for  flaw-seeded  specimens  of  7075-T73 
Forging  and  Ti-6A1-4V  beta  STOA  Forging  to  components  with  flaws  of  various  sizes  and 
shapes.  It  was  also  demonstrated  that  the  effective  stress  computed  using  the  notch  sensitivity 
factor  q  as  used  by  the  classical  methods  to  account  for  the  effects  of  notches  on  the  predicted 
fatigue  life  of  structural  components,  is  equivalent  to  computing  the  average  stress  over  a 
solution-dependent  domain  defined  by  J3.  This  approach  allows  the  extension  of  the  classical 
methods  developed  for  notches  to  flaws  of  arbitrary  shape. 

During  the  Option  3  (Task  6)  we  described  the  step-by-step  procedures  to  obtain  the  reference  S- 
N  curve  from  unnotched  NACA  specimens  of  24S-T3,  75S-T6  and  SAE  4130;  the  calibration 
procedures  to  obtain  the  J3-HSV  curves  for  each  material  utilizing  the  fatigue  data  from  notched 
NACA  specimens,  and  the  prediction  of  fatigue  life  based  on  the  calibration  curve,  including  the 
measure  of  the  difference  between  model  prediction  and  the  experimental  results. 

The  SimGov  Data  Management  Portal  (https://ftd.esrd.com)  was  updated  with  the  material 
developed  during  the  Task  6,  Option  3  Period,  as  shown  in  Appendix  4. 
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Appendix  1:  The  statistical  model  S2 


The  statistical  model  S2  differs  from  the  statistical  model  Si  only  in  that  the  standard  deviation  is 
assumed  to  be  a  function  of  aeq.  We  will  use  the  functional  form  proposed  by  Pascual  and 
Meeker  in  [3]: 

£  _  2QB1+B2  l°glO  ^eq  A'l 

Therefore  the  statistical  model  S2  is  characterized  by  six  parameters  02  =  Mi  A2  B\  B2  ju/Sf}.  The 
computed  parameter  values  are  shown  in  Figure  40. 


Figure  40:  Unnotched  NACA  data.  Computed  and  input  parameters  for  statistical  model  2.  Parameters  for 

the  design  curve  are  for  p =0.001. 

The  function  called  MMPDSmodel2RFL( mat ). m  is  the  main  routine  for  Model  S2.  Its  operation 
is  analogous  to  that  of  Model  Si.  The  PDF  and  the  CDF  functions  are  called  custpdfRFLI.m  and 
custcdfRFL2.m  respectively. 

For  Model  S2  the  seed  values  B\  and  B2  are  estimated  from 


Statistical  Model  S2 
NACA  materials:  RFL,  variable  STD 


24S-T3 

75S-T6 

SAE  4130 

A1 

7.344 

6.279 

7.177 

A2 

2.053 

1.484 

1.717 

B1 

-9.509 

1.137 

-6.083 

B2 

5.237 

-1.242 

2.946 

mu_f 

1.329 

1.430 

1.734 

s f 

0.0619 

0.0398 

f  0.0397 

Input  and  output  parameters 


24S-T3 

75S-T6 

SAE  4130 

U 

35 

35 

75 

L 

15 

15 

45 

s f 

0.05 

0.05 

0.05 

LL 

-572.85 

-879.89 

-477.30 

AIC 

1155.7 

1769.8 

964.6 

Parameters  for  the  design  curves 


24S-T3 

75S-T6 

SAE  4130 

A1 

14.458 

5.320 

5.261 

A2 

6.423 

0.972 

0.814 

A3 

-0.367 

21.200 

40.387 

log10  s  =  Bl+B2  log10  creq 


Letting 


log„A';-(A<'>-A<'»log,o(<7,s-A<'>)) 


And  Xj  =  log10  MM  =  log10  Sj  we  find  B\,  B2  by  linear  regression  using 

y j  ~  Bx  +  B2Xj,  j  =  l,2,...,n, 

The  computations  are  performed  using  the  Matlab  function  SeedValuesS2.m.  The  seed  values  for 
Model  S2  are: 

<&={#"  AS'>B,B2  log, 0(<>)  0.05} 


Page  5 1 


ESRD,  Inc. 


SBIR  II  Option  3  -  Task  6 


N68335-1 6-C-0088 


Discussion 

Since  Model  Si  is  a  special  case  of  Model  S2,  we  must  have  LL(  )  <  LL(§2).  This  is  found  to 
hold  numerically  for  the  three  materials  investigated.  The  AIC  values  take  into  consideration  the 
number  of  parameters  used  in  the  statistical  model.  The  AIC  values  are  consistently  smaller  for 
Si  than  for  S2,  therefore  S2  is  more  efficient  than  Si. 

Anomalous  behavior  is  exhibited  by  S2  for  the  24S-T3  and  SAE  4130  materials,  however.  One 
expects  the  standard  deviation  to  increase  with  decreasing  aeq  values  therefore  the  parameter  Bo 
in  Eq.  (A-l)  should  be  negative.  On  the  contrary,  we  find  that  for  these  two  materials  Bo  is 
positive  and  B\  is  negative  (see  Figure  40). 

It  was  possible  to  find  parameters  §2  such  that  Bo  was  negative  through  experimenting  with 
various  seed  parameter  values,  however  in  those  cases  LL(  Qx )  >  LL(§2)  which  clearly  indicates 
that  those  parameters  do  not  correspond  to  a  local  minimum.  It  is  also  surprising  that  the 
parameter  A3  of  the  design  curve  is  negative  for  24S-T3.  The  design  curve  shown  in  Figure  41 
for  Model  S2  is  comparable  with  the  design  curve  shown  in  Figure  1  for  Model  Si. 
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Figure  41:  S-N  data  for  24S-T3  aluminum  alloy.  The  mean  corresponding  to  Model  Si  is  labeled  P2(tfeq)’  the 

theoretical  design  curve  corresponding  top  =  0.001  is  labeled  d2(aeq).  Runouts  are  represented  by  the 

triangular  markers. 
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Appendix  2:  Comments  on  the  Random  Fatigue  Limit  Model 

On  page  279  in  reference  [4]  the  marginal  cumulative  distribution  for  the  normal-normal  model 
is  defined  as  follows: 


/X 

- $w\v 

-oo  ^7 


w  —  /.i(x,  v ,  9) 


a 


<f>v 


V  —  H-y 


'1 


dv 


A-2 


Where 


$ 


w\v 


w  —  n(x,  v,  9) 


w  ~  (/?o  +  ffi  log  [exp  (a;)  -  exp(t;)]’ 

o\pi 


A- 3 


and 


<pv 


A-4 


Let  us  consider  the  normal-normal  model.  From  Table  1  in  [4]  we  get: 

fio  =  30.272,  /3\  =  -5.1000,  a  =  0.289,  JuY  =  5.366,  ar  =  0.031 

Using  these  data  we  compute  the  marginal  CDF  for  various  stress  levels.  The  results  are  shown 
in  Figure  42.  In  the  neighborhood  of  s  =  exp  (/;,,)  «  214  MPa,  the  CDF  values  do  not  appear  to 
converge  to  1  as  w  ->oo.  The  reason  for  this  is  that  for  large  w  the  function  <E>wiv  is  unity 
independent  of  x  (where,  by  definition,  x  =  log  s).  On  the  other  hand,  the  integral  of  <p\j  depends 
on  x.  Therefore  for  w->oowe  have 


F{w,  x,  9)  —¥  g{x,  9)  :=  J  — exp  ^  dv.  A-5 

The  computed  values  of  g(x,  6)  are  shown  in  Table  1  and  displayed  in  Figure  42.  It  is  seen  that 
g(x,  0)  is  sensitive  to  s  when  x  is  close  to  jUr 

In  conclusion,  the  marginal  CDF  and  PDF  functions  should  be  scaled  by  l/g(x,  6)  to  ensure  that 
the  CDF  is  unity  for  large  w,  independent  of  s.  This  scaling  affects  the  definition  of  the 
likelihood  function  and  hence  6. 
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Region  of  interest 


0.8136 


0.5000 
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Figure  42:  Cumulative  distribution  functions  corresponding  to  the  indicated  stress  levels. 


Table  12:  The  computed  values  of  g(x,  0)  corresponding  to  the  stress  values  shown  in  Figure  42.  The  stress 

values  are  in  MPa  units. 


s  =  350  s  =  300  s  =  250  .s  =  220  s  =  214  s  =  210 

g(x,0) 

1.0000  1.0000  1.0000  0.8136  0.5000  0.2711 

Correspondence  on  this  point  with  Professors  Pascual  and  Meeker  is  reproduced  in  the 
following. 

Correspondence  with  Professors  F.  G.  Pascual  and  W.  Q.  Meeker 

Date:  Sat,  12  Nov  2016  13:15:23  +0000 
From:  "Szabo,  Barna"  <szabo@wustl.edu> 

To:  "wqmeeker@iastate.edu"  <wqmeeker@iastate.edu>, 

"jpascual@math.wsu.edu"  <jpascual@math.wsu.edu> 

Subject:  Question  relating  to  the  random  fatigue  limit  model 

Dear  Professors  Pascual  and  Meeker, 

With  reference  to  your  joint  paper  published  in  Technometrics  in  1999,  I  believe  that  the 
marginal  PDF  and  CDF  functions  defined  in  that  paper  should  be  scaled  as  explained  in  the 
appended  note. 

I  will  greatly  appreciate  if  you  would  let  me  know  whether  you  agree  with  the  conclusion. 
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With  best  regards, 

Bama  Szabo 

From:  Jave  Pascual  <jpascual@math.wsu.edu> 

Sent:  Monday,  November  14,  2016  1:10  PM 
To:  Szabo,  Barna 

Cc:  Jave  Pascual;  wqmeeker@iastate.edu 

Subject:  Re:  Question  relating  to  the  random  fatigue  limit  model 

Professor  Szabo, 

Thanks  for  interest  in  the  RFL  model.  Below  are  some  thoughts  related  to  your  inquiry: 

One  of  the  motivations  for  the  RFL  model  was  to  "model"  infinite  life  whenever  the  unit’s 
fatigue  limit  exceeded  applied  stress.  One  way  of  looking  at  it  is  that  there  is  a  positive 
probability  mass  at  infinite  life  and  that  probability  is  the  same  as  the  probability  that  a  unit’s 
fatigue  limit  falls  above  the  applied  stress.  Infinite  life  and  the  "probability  deficit"  are  still 
accounted  for  in  the  likelihood  when  we  compute  survival=  1-cdf. 

Also,  when  you  normalize  or  standardize  the  pdf  or  cdf,  you  convert  marginal  to  conditional 
probabilities  wherein  you  are  making  the  additional  assumption  that  the  fatigue  limit  is  lower 
than  the  applied  stress  FOR  ALL  units. 

The  feature  of  placing  a  distribution  on  the  fatigue  limit  allows  the  model  to  mechanistically  (as 
opposed  to  curve-fitting)  describe  the  curvature  that  we  often  see  in  the  S/N  curve. 

Jave 

F.  G.  Pascual 
Associate  Professor 

Department  of  Mathematics  and  Statistics 
Washington  State  University 
PHONE:  (509)335-3126 
FAX:  (509)335-8369 
E-MAIL:  iave@wsu.edu 

Meeker,  William  Q  [STAT]  <wqmeeker@iastate.edu> 

Dear  Professor  Szabo: 

The  distribution,  as  written  in  the  paper  is  correct.  A  feature  of  the  RFL  model  is  that  there  is 
always  a  positive  probability  that  life  will  extend  to  infinity.  Think  of  it  this  way.  The  fatigue 
limit  is  a  random  variable,  but  once  it  is  realized  on  a  given  unit,  it  is  fixed.  If  the  applied 
constant  stress  for  that  unit  is  less  than  the  fatigue  limit,  then  a  crack  will  never  initiate. 
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Best  regards, 

Bill  Meeker 

Remarks 

If  the  conditional  PDF  is  not  scaled  then,  as  seen  in  Figure  42,  when  the  applied  stress  is  equal  to 
ju7  (214  MPa)  then  there  is  a  50%  chance  for  infinite  life.  This  does  not  seem  plausible. 

Scaling  the  conditional  PDF  means  that  failure  will  always  occur.  Using  the  scaled  PDF  should 
be  considered  as  an  alternative  to  the  RFL  model  proposed  by  Pascual  and  Meeker  in  [4], 
Pascual  proposed  yet  another  variant  of  the  RFL  model  in  [7], 

As  noted  by  Professor  Pascual,  scaling  converts  marginal  probability  to  conditional  probability 
implying  the  assumption  that  the  fatigue  limit  may  be  zero.  As  a  practical  matter  this  is  not  a 
concern  because  the  calibration  data  would  not  justify  extrapolation  beyond  100  million  cycles. 

We  should  note  also  that  aluminum  alloys,  unlike  ferrous  and  titanium  alloys,  do  not  have  a 
fatigue  limit.  Therefore  the  fatigue  limit  model,  applied  to  these  materials,  should  not  be 
considered  valid  for  fatigue  cycles  exceeding  some  limit  which  is  not  greater  than  108 
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Appendix  3:  Configuration  of  the  specimens  identified  in  Table  2 

NACA  specimen  dimensions  and  corresponding  finite  element  meshes  used  for  the  analysis. 
•  Specimen  #  1  -  Unnotched  dogbone  specimen  -  Actual  Kt  =  1.00 


•  Specimen  #2  -  Hole  type  notch  -  Actual  Kt  =  2.11 


mmw 


•  Specimen  #3  -  Edge-cut  notch  -  Actual  K,  =  2. 17 
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•  Specimen  #4  -  Fillet-type  notch  -  Actual  Kt  =  2. 19 
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Specimen  #5  -  Edge-cut  notch  -  Actual  K,  =  4.43 
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•  Specimen  #6  -  Fillet-type  notch  -  Actual  K,  =  4.83 
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•  Specimen  #8  -  Edge-cut  notch  -  Actual  Kt=  1 .62 


utmu 

HU  WU 


•  Specimen  #9  -  Edge-cut  notch  -  Actual  K,  =  4.48 
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[f 


mm 


R  =  0.0035  ±  0.0005 


•  Specimen  #10  -  Edge-cut  notch  -  Actual  Kt  =  4.41 
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Appendix  4:  Updates  to  SimGov  Portal  during  the  Task  6  of  the 
Option  3  Period 

All  files  for  Task  6  referenced  throughout  this  report  are  located  in  the  SimGov  Data 

Management  Portal  https.V/ftd. esrd. com,  under  the  Data  Group  NACA  Data.  The  file  names  and 

contents  are  as  follows: 

1.  Model  SI  Archived.zip  and  Model  S2  Archived.zip  -  Matlab  functions  and  input  data  for 
statistical  models  Si,  and  Si,  according  to  the  details  given  in  Section  2.3  and  in  Appendix  1. 

2.  naca-tn-24S-T3-(Raw  Datal.xlsx  -  Fatigue  data  of  93  unnotched  specimen  extracted  from 
NACA-TN-2324,  Table  2  is  tabulated  in  this  Excel  file,  worksheet  “Kt=l  data”.  Of  the  93 
specimens,  10  were  runouts  and  11  were  disqualified  tests. 

3.  naca-tn-75S-T6-(Raw  Datal.xlsx  -  Fatigue  data  of  86  unnotched  specimen  extracted  from 
NACA-TN-2324,  Table  3  is  tabulated  in  this  Excel  file,  worksheet  “Kt=l  data”.  Of  the  86 
specimens,  12  were  runouts  and  3  were  disqualified  tests. 

4.  naca-tn-SAE-4130-(Raw  Datal.xlsx  -  Fatigue  data  of  59  unnotched  specimen  extracted  from 
NACA-TN-2324,  Table  4  is  tabulated  in  the  Excel  file,  worksheet  “Kt=l  data”.  Of  the  59 
specimens,  16  were  runouts  and  6  were  disqualified  tests. 

5.  BetaSummary24S-T3.xlsx  -  Values  of  the  parameter  (3  computed  for  all  test  cases  of  the 
eight  specimen  types  used  for  calibration,  and  for  values  of  a  =  0.0,  0.2,  0.4.  0.5,  0.6,  0.8,  1.0. 
Material  aluminum  24S-T3. 

6.  BetaSummary75S-T6.xlsx  -  Values  of  the  parameter  f3  computed  for  all  test  cases  of  the 
eight  specimen  types  used  for  calibration,  and  for  values  of  a  =  0.0,  0.2,  0.4,  0.5,  0.6,  0.8,  1.0. 
Material  aluminum  75S-T6. 

7.  BetaSummarySAE-4130.xlsx  -  Values  of  the  parameter  (3  computed  for  all  test  cases  of  the 
eight  specimen  types  used  for  calibration,  and  for  values  of  a  =  0.0,  0.2,  0.4,  0.5,  0.6,  0.8,  1.0. 
Material  steel  SAE-4130. 

8.  FittedBetaSummary.xlsx  -  Summary  of  /7-values  computed  from  Eq.  (30)  and  Eq.  (31)  for 
specimen  9,  all  materials  (see  Table  3);  summary  of  Bayes  factors  for  specimen  9,  all 
materials  (see  Table  4);  summary  of  optimal  values  of  a  for  each  material  and  /7-fit  (see 
Table  11).  Macro  Excel  files  used  for  the  computations  of  the  (3  values  for  each  test  case  of 
the  8  specimens  types  used  for  calibration  are  in  the  following  files  (with  xx  =  00,  02,  04,  05, 

06,  08,  10): 

o  Prediction  (24S-T3  SWT)-RFL-Model- 1  -Alpha-xx.xlsm 

o  Prediction  (75S-T6  SWT)-RFL-Model-l -Alpha-xx.xlsm 

o  Prediction-S  AE-4 1 30-RFL-Model- 1  -  Alpha-xx.xlsm 
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9.  RankingData-24S-T3-Specimen-9.xlsx  -  Ranking  of  /?  prediction  using  Eq.  (26)  to  compute 
G^ ,  given  ~p(arHSA)  and  j3 AVG (a j,HS A),  for  each  of  the  test  records  corresponding  to 
specimen  #9.  Material  aluminum  24S-T3. 

10.  RankingData-75S-T6-Specimen-9.xlsx  -  Ranking  of  f3  prediction  using  Eq.  (26)  to  compute 
G(^ ,  given  j3(cCj,HSA)  and  f3AVG(cCj,HSA), for  each  of  the  test  records  corresponding  to 
specimen  #9.  Material  aluminum  75S-T6. 

11.  RankingData-SAE-4130-Specimen-9.xlsx  -  Ranking  of  /?  prediction  using  Eq.  (26)  to 
compute  G^ ,  given  P(cCj,HSA)  and  j3 AVG (a j,HS A),  for  each  of  the  test  records 

corresponding  to  specimen  #9.  Material  steel  SAE-4130. 

12.  Calibration-All-(24S-T3-SWT)-RFL-Model-1-Alpha-O2.xlsm  and  Calibration- A11-(24S-T3- 
S WT)-RFL-Model- 1  -  Alpha-08  .xlsm  -  Updated  J3-HSA  calibration  curves  for  all  nine 
specimen  types  for  the  optimal  values  of  a.  Material  aluminum  24S-T3. 

13.  Calibration-All-(75S-T6-SWT)-RFL-Model-1-Alpha-O6.xlsm  and  Calibration- A11-(75S-T6- 
S WT)-RFL-Model- 1  -  Alpha- 1 0.xlsm  -  Updated  J3-HSA  calibration  curves  for  all  nine 
specimen  types  for  the  optimal  values  of  a.  Material  aluminum  75S-T6. 

14.  Calibration- All-SAE-4 1 30-RFL-Model- 1  - Alpha-OO.xlsm  -  Updated  J3-HSA  calibration 
curves  for  all  nine  specimen  types  for  the  optimal  value  of  a.  Material  steel  SAE-4130. 

15.  BetaZ-Score- Alpha-0. 8(24S-T3).xlsx  -  Detailed  computation  of  the  frequency  plots 

described  in  Section  2.5.1.  Material  aluminum  24S-T3  for  a  =  0.8. 

16.  BetaZ-Score- Alpha-1.0(75S-T6).xlsx  -  Detailed  computation  of  the  frequency  plots 

described  in  Section  2.5.1.  Material  aluminum754S-T6  for  a-  1.0. 

17.  BetaZ-Score- Alpha-0.0(SAE-4130).xlsx  -  Detailed  computation  of  the  frequency  plots 
described  in  Section  2.5.1.  Material  steel  SAE-4130  for  a  =  0.0. 
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